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SECTION  I 


INTRODUCTION 

The  optical  properties  of  a  medium  can  be  altered  by  a  laser  beam  of 
sufficient  intensity.  This  alteration  can  range  from  relatively  small  changes 
in  the  index  of  refraction  leading  to  diffraction  effects  to  complete  breakdown 
of  the  propagating  medium.  Thus,  the  interaction  between  a  laser  beam  and 
the  surrounding  environment  can  impose  limitations  on  the  size  and  intensity 
of  the  beam  that  is  to  be  used  for  a  given  purpose. 

In  the  case  of  air,  laboratory  experiments  with  the  10.  6  n  radiation 

from  a  CO  laser  have  shown  that  atmospheric  breakdown  occurs  at  an  inten- 
2  9  2 

sity  of  approximately  2  X  10  watts/cm  for  sea  level  density  (Ref.  1).  It  is 
shown  in  Section  II  that  thi3  is  about  the  same  intensity  predicted  by  theoret¬ 
ical  considerations  if  it  is  assumed  that  initially  there  is  one  free  electron 
3 

per  cm  . 

Further  experiments  indicated  that  the  breakdown  threshold  for  air 
8  2 

could  be  as  low  as  10  watts/cm  (Ref.  2).  Since  the  later  experiments  were 
carried  out  with  air  under  normal  conditions  while  the  earlier  experiments 
used  cleaner  air  meeting  laboratory  standards,  it  was  postulated  that  the 
lowering  of  the  breakdown  threshold  was  due  to  the  presence  of  an  increased 
amount  of  contaminants  in  the  dirty  air. 

Contaminants  such  as  quartz  (SiO?)  and  sapphire  (Al.OJ  are  assumed 
to  be  present  in  the  form  of  solid  particles  ranging  in  size  from  1  ^  to  100  jj 
radius.  These  particles  absorb  energy,  finally  reaching  the  vaporization 
temperature  and  moving  out  into  the  air.  Because  of  the  higher  absorption 
of  energy  the  contaminants  can  provide  the  initial  source  of  electrons  lead¬ 
ing  to  atmospheric  breakdown. 

The  purpose  of  the  present  study  is  to  develop  a  model  suitable  for 


investigation  of  the  effects  of  contaminants  on  atmospheric  breakdown.  This 
model  utilizes  where  possible,  the  known  material  properties  of  the  contam¬ 
inant  particle  to  determine  the  absorption  of  the  laser  energy  and  the  subse¬ 
quent  heating.  A  one  dimensional  Lagrangian  hydrodynamic  code,  CELaB 
(Contaminant  Effects  on  Laser  Air  Breakdown),  has  been  developed  to  com¬ 
pute  the  motion  of  the  vaporized  material. 

Semi-emperical  expressions  for  the  relevant  molecular  cross  sections 
are  developed  to  describe  the  interactions  between  the  free  electrons  from 
the  contaminant  particle  and  the  atmosphere.  The  calculated  values  for  the 
cross  sections  are  compared  with  the  latest  data  available  and  found  to  give 
reasonable  results  over  a  wide  range  of  electron  energy. 

The  results  of  preliminary  calculations  on  quartz  (SiO  )  and  sapphire 
(•^2^3^  particles  of  varying  sizes  indicate  that  atmospheric  contaminants  are 
an  important  factor  in  the  laser  induced  breakdown  of  air.  First,  they  supply 
a  substantial  density  of  free  electrons  for  the  coupling  of  laser  energy. 
Secondly,  because  of  the  relatively  high  density  of  the  solid  contaminant,  the 
local  degree  of  ionization  in  the  neighborhood  of  the  evaporated  particle  will 
be  high  making  the  electron-ion  interaction  the  dominant  elastic  collision 
process.  This,  in  turn,  increases  the  rate  at  which  an  electron  gains  energy 
from  the  laser  and  decreases  the  breakdown  threshold  intensity. 
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SECTION  II 


AIR  BREAKDOWN  ANALYSIS 

The  intensity  of  high-power  laser  beams  propagating  through  the 

atmosphere  is  limited  by  atmospheric  breakdown.  For  the  10.6-micron 

radiation  of  the  CO.  laser,  experiments  on  ordinary  laboratory  air  show 

9  2 

that  breakdown  is  initiated  at  laser  powers  of  about  2  x  10  watts/cm  .  If 
it  is  assumed  that  there  is  an  initial  density  of  free  electrons  in  the  focal 
volume  of  the  laser  beam  then  the  breakdown  process  can  be  described  by 
electron  cascade  theory  (Ref.  3).  In  this  theory,  the  free  elections  in  the 
focal  volume  of  the  beam  absorb  energy  by  inverse  bremsstrahlung.  When 
the  rate  of  energy  absorption  is  large  enough  for  an  electron  to  reach  the 
ionization  energy  during  the  laser  pulse  duration,  there  follows  a  net  rate 
of  growth  for  the  free  electron  population.  If  this  growth  rate  exceeds  the 
rate  of  loss  of  electrons  to  attachment,  diffusion,  etc. ,  then  the  electron 
density  increases  with  increasing  time  and  breakdown  of  the  gas  is  taking 
place.  The  breakdown  of  the  gas  is  complete  when  the  growth  rate  is  large 
enough  that  complete  ionization  takes  place  within  the  time  of  the  pulse  dur¬ 
ation.  The  process  can  be  described  by  the  rate  equation, 

t 

-  aL)dt 

Ne(t)  =  Ne(0)e  (1) 

where  Ng(t)  is  the  time  dependent  electron  density,  y  is  the  growth  rate, 
and  is  the  loss  rate.  Assuming  that  the  rates  are  constant  over  the  dur¬ 
ation  of  the  laser  pulse,  the  criterion  for  breakdown  requires  that  the  growth 
rate  satisfies  the  condition 

7  >«L  (2) 
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At  standard  temperature  and  pressure  conditions  (STP)  in  air,  complete 
breakdown  requires 

45 

v’a L*r  <3> 

p 

where  it  has  been  assumed  that  a  nominal  value  of  N  (0)  =  1  cm  exists 

19  -3  e 

and  N^ftp)  =  3  X  10  cm  In  the  remainder  of  the  report,  the  term 
"breakdown"  will  always  refer  to  complete  breakdown  unless  otherwise 
specified. 

The  dominant  loss  mechanism  in  air  is  3-body  attachment  to 
molecular  oxygen, 


02+M+e«S02+M 


(4) 


where  M  is  a  third  body.  The  rate  for  this  process  is  (Ref.  4) 


~  7.  5  X  107  (p/p  sec  1 


(5) 


where  pQ  is  the  sea  level  density.  At  STP  conditions, 


-  7  45  -1 

y  ~  7.  5  x  lo'  +p.  8ec 

p 


(6) 


Thus,  for  very  long  pulses,  the  threshold  growth  rate  for  breakdown  i6 


y  =  7.  5  xlO7  sec-1 


(7) 


while  for  pulses  in  the  nanosecond  range 


-  45  -1 

y  sec 

P 


(8) 
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.  Schultz 


Electron  Energy  -  ev 


Figure  1.  Sum  of  ^  Vibrational 
Cron  Sections 


1.  CROSS  SECTIONS 


a.  Vibrational  Excitation 

The  vibrational  excitation  of  N2  by  electron  impact  has  been 
studied  quite  extensively  both  experimentally  and  theoretical!/  (Ref.  5). 

In  the  theoretical  analysis,  it  has  been  assumed  that  the  vibrational  excita¬ 
tion  takes  place  through  the  formation  of  a  compound  negative  ion  state. 
Figure  1  shows  the  total  cross  section  for  the  excitation  of  the  first  eight 
vibrational  states.  It  should  be  noted  that  the  interaction,  having  a  peak  at 
an  energy  of  ~  2.  3  ev,  is  a  resonant  interaction  having  a  cut  off  at  ~3.  5  ev. 
For  the  calculations  to  be  carried  out  in  this  report,  an  adequate  represen¬ 
tation  can  be  obtained  by  considering  the  cross  section  as  a  constant 
av  ~  3  *  10  over  the  energy  range  1.9  £  f  £  2.  8  ev. 

b.  Electronic  Excitation 

The  determination  of  the  absolute  value  of  the  cross  sections  for 
the  excitation  of  molecular  electronic  levels  poses  a  difficult  problem  both 
experimentally  and  theoretically.  The  laost  frequently  quoted  values  for 
N2  and  are  those  given  by  Engelhardt,  Phelps,  and  Risk  (Ref.  6)  for  N2 
and  Hake  and  Phelps  (Ref.  7)  for  (X,.  The  method  used  by  these  authors  is 
to  numerically  solve  the  time-and  space -independent  Boltzmann  equation 
for  the  electron  energy  distribution  function.  Averages  over  this  distribu¬ 
tion  function  determine  transport  coefficients  which  are,  in  turn,  compared 
with  experiment.  Appropriate  adjustments  are  made  to  the  cross  sections 
until  a  satisfactory  fit  with  the  experimental  data  for  the  transport  coeffi¬ 
cients  are  obtained.  The  thresholds  for  the  various  cross  sections  are 
obtained  from  other  experimental  data.  The  cross  sections  giving  the  best 

fit  to  the  experimental  transport  coefficients  over  a  wide  energy  range  are 
then  plotted. 

More  recently,  Borst  (Ref.  8)  has  made  an  experimental  deter¬ 
mination  of  the  absolute  cross  sections  for  the  excitation  of  the  A^E+, 


1  3  + 

a  II  and  £  2  states  of  N-«  Ajello  (Ref.  9)  has  also  made  measurements 

•  •  6  j 

of  the  cross  section  for  the  excitation  of  the  a  IK  state  of  N,.  There 

8 

appears  to  be  substantial  disagreement  between  the  two  measurements.  The 
experimental  data  will  be  shown  in  appropriate  portions  of  this  section 
where  comparisons  are  to  be  made  with  calculated  values  of  the  cross 
sections. 

An  approximate  method  for  the  calculation  of  the  relevant  cross 
sections  can  be  obtained  from  classical  mechanics.  The  transfer  of  energy 
between  two  electrons  is 


\ 


Ac  =  C  «in  6/2 


where  c  is  the  kinetic  energy  of  the  incident  electron  and  6  is  the  scatter- 
ing  angle.  The  relation  between  the  impact  parameter ,  energy,  and  scat¬ 
tering  angle  is  (Ref.  10) 

2 

b=~-cote/2  (10) 


hence. 


Ac  = 


1  + 


(H) 


where  b  is  the  impact  parameter  and  e  is  the  electronic  charge. 

According  to  classical  mechanics,  the  cross  section  is  defined 
by  the  integral 


a  (c)  =  2ir 


M 

/ 


bdb  =  ff  b 


M 


(12) 


where  b^  is  the  maximum  value  of  the  impact  parameter  for  the  particular 
interaction  process  under  consideration.  For  electronic  excitation,  the 


7 


required  value  of  the  impact  parameter  can  be  determined  from  equation 
(11)  by  setting  Ac  =  ^  . 


f 


—  «  -  *,) 


(13) 


and 


**«> 


ffe 


(c  -  cx) 


(14) 


where  c  is  the  excitation  energy.  Equation  (14)  has  a  peak  value  at  c  =  2c 

X  3 

giving  a  peak  value  of  the  cross  section 


a  (2*  )  = 

X  *x 


ffe 


1 6c 


(15) 


Equation  (14)  seems  to  give  reasonable  behavior  for  singlet- 
singlet  transitions  such  as  -*  a*n  .  However,  for  the  singlet-triplet 

If  3  ^ 

transitions  such  as  Z  ■♦A  equation  (14)  does  not  fall  off  rapidly  enough 

o 

with  increasing  energy.  The  functionaJ><form  for  the  triplet  excitation  cross 
section  is 


0x(c>  ^OjCatc  -  cx)}e 


-ac 


(16) 


where  Oj  and  a  are  constants  to  be  determined. 

The  maximum  value  for  the  cross  section  given  by  equation  (16) 
occurs  at  the  energy 


’  7  +  fx 
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Setting  =  2f^  so  that  the  peak  cross  section  occurs  at  the  same  energy 
as  predicted  by  equation  (14)  gives  the  result  a  *  cx  .  and 


<*x<c)  =  oL 


(17) 


Setting  (:2(  in  equation  (17)  and  equating  the  resulting  cross  section  to 
the  peak  cross  section  given  by  equation  (14)  leads  to  an  expression  for  . 


°1 


0.461 


(18) 


and 


(19) 


Equations  (14)  and  (19)  will  be  used  to  calculate  the  various  electronic  exci- 
tation  processes  of  interest  in  this  report  and  the  results  compared  with 
existing  experimental  data  and  the  results  of  other  calculations. 

3  + 

(1)  Nitrogen  Cross  Section  for  A  Zu  State 

The  excitation  energy  for  this  state  is  about  Cx  =  8.  7  ev 

(Ref.  11);  hence,  the  calculated  peak  value  of  the  cross  section  will  occur 

at  (  =  13.4  ev  and  the  peak  value  of  the  cross  section  calculated  from 

- 17  2 

equation  (19)  is  0^(2^)  =  9*  1  X  10  cm  . 

The  experimental  data  of  Borst  (Ref.  8)  place  the  peak 
at  c  ~  11  ev,  which  is  in  reasonable  agreement  with  calculated  value.  The 
measured  value  of  the  peak  cross  section  is  <jx  s  ^5.  3*  ^  j  X  10  ^cm^. 

The  upper  limit  is  very  close  to  the  calculated  value. 
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The  calculated  peak  value  of  Engelhardt,  et  al.  (Kef.  6)  is 
- 1 6  2 

0^  M  0.  6  x  10  cm  ,  but  their  peak  occurs  at  about  7.  5  ev,  whicu  is  about 

4  ev  lower  than  the  experimental  value.  The  calculated  value  of  Chung  and 

and  Lin  (Ref.  12)  has  a  peak  value  at  about  11  ev,  but  their  calculated  peak 

-16  2 

cross  section  is  a  w  1  21  x  10  cm  . 
x 

The  cross  section  used  by  Kroll  and  Watson  (Ref.  11)  for 
3  + 

excitation  of  the  N,(A  £  )  level  is 

2  u 

CT^(e)  =  3.45  x  10  ^  (c  -  6.  7)  e  c.a^  (20) 

where  f  is  in  units  of  ev.  This  is  to  be  compared  with  the  equation  used  in 
this  report, 

•  ox(c)  =  1.01  X  lo"16  (C  -  6.7)  e‘  c/6'7  cm2  (21) 

Thus,  Kroll  and  Watson  obtain  a  peak  cross  section  about  a  factor  of 
3  higher  than  obtained  in  the. present  calculation.  This  is  attributed  to  the 
fact  that  Kroll  and  Watson  multiplied  by  a  scale  factor  to  give  the  measured 
low-frequency  breakdown  threshold  in  air. 

3  + 

Various  excitation  functions  for  the  A  £  state  of  N_  are  com- 

u  2 

pared  in  figure  2.  The  peak  value  error  bar  given  by  Borst  is  shown  to 
indicate  the  magnitude  of  the  experimental  uncertainties  involved.  The 
values  of  Kroll  and  Watson  equation  (20)  are  not  plotted  since  they  are 
similar  to  those  for  equation  (21)  except  off  scale  a  factor  of  about  3  higher. 

3 

(2)  Nitrogen  Cross  Section  for  B  State 

The  excitation  energy  for  this  state  is  about  7.4  ev,  giv- 

-17  2 

ing  a  peak  value  of  the  cross  section  at  14.  8  ev  of  7.  5  x  10  cm  calcu¬ 
lated  from  equation  (19).  Chung  and  Lin  (Ref.  12)  report  both  a  corrected 
experimental  excitation  function  for  and  a  theoretical  one.  These  are 

r  * 

.10 


v 


/  \ 

I  \ 


Engelhardt,  Phelps,  and 
Risk 


Chung  and  Lin 


Figure  2.  Excitation  Functions  for  the 
A3!  *  State  For  Ng 


shewn  in  figure  3.  Also  shown  is  the  theoretical  curve  obtained  from  equa¬ 
tion  (19)  which  is  numerically 


ax(c>  =  7.  16  X  10-l7(f  -  7.4)  e‘  c/7,4 


(3)  Nitrogen  Cross  Section  for  C3n  State 
- - - - u - 

The  excitation  energy  for  the  C  IIu  state  is  about  11.  2  ev 
giving  a  peak  value  of  the  cross  section  at  22.  4  ev  of  3.  3  x  I0_l7ev.  Chung 
and  Lin  (Ref.  12)  report  excitation  functions  fpr  this  state  both  experimen¬ 
tal  and  theoretical  from  other  sources  as  well  as  their  own  theoretical 
values.  These  are  shown  in  figure  4  along  with  the  experimental  values 
reported  by  Engelhardt,  Phelps  and  Risk  (Ref.  6).  For  this  case,  equation 
(19)  becomes  numerically 

ox(f)  =  2.  15  x  10"17(€  -  11.2)  e"  c/11,2  (23) 

values  from  which  are  shown  for  comparison  in  figure  4.  It  can  be  seen 
that  in  this  instance  the  peak  appears  somewhat  sharper  at  lower  energies 
than  given  by  the  classical  form  of  equation  (23). 

(4)  Nitrogen  Cross  Sections  for  the  a*n  State 

g - 

The  excitation  energy  for  the  a^  state  is  about  8.4  ev. 

g 

For  this  state,  as  discussed  above,  the  form  of  equation  (14)  is  appropri¬ 
ate.  The  peak  cross  section  in  this  instance  occurs  at  16.  8  ev  with  a  value 
-17  2 

of  5.  8  x  10  cm  .  Chung  and  Lin  report  a  theoretical  excitation  function 
for  this  state,  and  experimental  values  are  reported  by  Engelhardt,  Phelps 
and  Risk  (Ref.  6),  by  Ajello  (Ref.  9)  and  by  Borst  (Ref.  8),  all  of  which  are 
shown  in  figure  5.  The  numerical  form  of  equation  (14)  for  this  case  is 


0  (c)  =  1.94  x  10’l5(c  -  8.4)/c2 
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——  Equation  (24) 

— —  -  -  Chung  and  Lin 

—  —  Engelhardt,  Phelps, 

and  Risk 


0  Ajello 


Values  obtained  from  this  equation  are  also  shown  in  figure  5  for  compari¬ 
sons. 


c.  Ionization 

Extensive  data  exists  for  ionization  cross  sections  in  both  nitro¬ 
gen  and  oxygen.  Values  obtained  from  a  summary  of  these  data  are  shown 
in  figures  6  and  7.  Approximate  empirical  forms  wnich  fit  these  data  rea¬ 
sonably  well  at  least  over  the  peak  are  for  N  : 

2 


a  (c)  =  1.76  x  10 


-15 


(25) 


where  c  =  15.  58  ev,  the  ionization  threshold  for  N  .  For  O,  the  empirical 
equation  is: 


a. 

i 


(f)  =  5.  10  x  10“16 


e/l5c. 


(26) 


where  c  .  =  12.  06  ev,  the  ionization  threshold  for  O  .  Values  obtained 

Cm 

from  these  two  empirical  equations  are  shown  for  comparison  in  figures 
6  and  7. 

d.  Total  Scattering 

The  total  N2  scattering  cross  section  reported  by  Golden  (Ref. 
13)  is  shown  in  figure  8.  For  comparison,  an  empirical  fit  is  also  shown 
calculated  from  the  formula. 

<.TU).[l+2«-2<‘-2-5>2],10-15  (27, 

All  of  the  empirical  formulae  described  above  are  approximate, 
but  considering  the  large  uncertainties  present  in  much  of  the  experimental 
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Fitted  -  Eq.  (25) 
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Figure  7.  Total  Ionization  Cross 
Section 
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Figure  8 


Total  Electron  Scattering 
Cross  Section 


data  available,  these  formulae  are  useful  for  purposes  of  estimation. 

2.  BREAKDOWN  ANALYSIS 

An  analysis  for  the  determination  of  the  threshold  field  leading  to 
breakdown  will  be  given  in  this  subsection.  The  method  for  solving  the 
kinetic  equation  is  similar  to  the  method  used  by  Baraff  and  Buchsbaum 
(Ref.  14). 

An  approximate  kinetic  equation  for  electrons  gaining  energy  at  a 
rate  k  is 


(*t  H  k  *7 ) 1  =  *<*>  -  vxif^*  (28) 

where  Vxi  ^  and  v.  are  the  excitation  and  ionization  collision  fre¬ 

quencies,  and 


(Vx  +  2vi)f(c.  t)dc 


(29) 


Ihe  factor  of  2  enters  in  equation  (29)  since  electrons  are  gained  in  the  dis¬ 
tribution  following  an  ionization  event  in  which  the  ionizing  electron  is  left 
with  energy  f  and  ionizing  events  in  which  the  ionized  electron  has  an  energy 
C  •  The  approximation  for  the  kinetic  equation  assumes  that  the  rate  at 
which  an  electron  gains  energy  in  the  oscillating  electric  field  can  be  repre¬ 
sented  in  terms  of  an  "effective"  dc  field.  This  approximation  has  been 

shown  to  be  reasonable  if  the  radiation  frequency  is  much  larger  than  electron 
growth  rate  (Ref.  14). 

Defining  the  Laplace  transforms  of  the  distribution  function  f(c,t)  and 
the  function  g(t). 
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• 

• 

He.  p)  =  [  f(t,t)  e"pt  dt 
•'o 

(30) 

c 

m 

g(p)  =  [  g(t)  e"pt  dt 
•'o 

• 

(31) 

s  v(c)?U,  p)df 
"o 

where  v(c)  ■  2v^(c)  +  vjc)*  Equation  (28)  becomes 

(k  £  +  p  +  vxi)  r<tlP>  *  «<<>  c 1  +  «  (pH 

(32) 

• 

where  the  assumption 

■* 

f(c.°)  =  fi(e) 

(33) 

has  been  used.  Dividing  equation  (32)  by  [1  +  g  (p)]  , 

(*  ^•+p  +  vxijef(6p)  =  j(c> 

(34) 

where 

Su.  P) . 

1  +  g(p) 

(35) 

Multiplying  equation  (35)  by  v(f)  and  integrating  over  all  t 
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m 

,(<,eu,p,<l‘*  T^T  [  «<>?<<•  e'*' 


Using  equation  (31),  equation  (36)  can  be  rewritten  as 


|(p)  =  -5i£l - 

i  -  fir  (pi 


where 


ff<P)  ■  /  v(c)ff(c.p)df 


Taking  the  Laplace  inversion  of  equation  (37), 


£  +  i 


J—  [  pt 

2,1  J  l-«T(p)  P 


fi‘  i- 


Equation  (39)  can  be  integrated  by  closing  the  contour  in  the  left  half  of  the 
complex  p- plane  giving 


g(t)  =  - 


V 71 

m \ 

ydp  ^p  = 


where  y  is  the  ~ue  of  p  for  which  ff(p  =  y)  =  1.  Thus,  y  is  the  growth  rate 
for  electron  p  ■  Auction  and  is  determined  by  the  condition  ff(y)  =  1. 

Solving  e,  nation  (34)  for  G((,  y), 
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S(t.y)  = 


( 

E«  ° 


(41) 


Substituting  equation  (41)  into  equation  (38)  and  using  the  relation  ff(y)  =  1  , 


C /  v(<>  “p  [-  £  (**  *f  *«*)]  d*  * 


Equation  (42)  can  be  written  in  the  equivalent  form 


dc 


•  c 

+  £  f  vi(t,exp  [-7  (»'*+/  vxiu'|d<')]  d< 
O  o 

•  ( 

-  k  /  “P  [-  £  (*«  *j  Vxl(«')dc')]d€  ■  1 


The  first  integral  of  equation  (43)  integrates  to  unity,  hence, 

.e 

j  / 


/  Vj(c)  « 
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k 


•  E  /  vxi(<‘ 


>d( 


_  J_o _ _ 

r'-fi  . 

I  •  Jo  df 


(42) 


(43) 


(44) 


When  the  energy  is  less  than  the  ionization  energy,  the  numerator 
of  equation  (44)  is  zero,  therefore  equation  (44)  can  be  rewritten  as 


-  v  Ac/k  • 

X  f 


b 


/ 

k  k  / 
U>  e  \ 


* 

£  J  vx i(»',d'' 


*0 


where 


Vx  ^  d€# 


and  €•»  are  the  ionization  and  excitation  energies,  Respectively. 

Equation  (45)  can  be  written  in  terms  of  a  mean  value  of  the  ioniza¬ 
tion  frequency,  <- 


y~  <  Vj>  e 


-  Vx  Ac/k 


[-*] 


where  B  is  a  correction  factor  describing  the  deviation  of  the  ionization  fre¬ 
quency  from  its  mean  value.  The  value  of  B  can  range  from  B  =  0;  corres¬ 
ponding  to  no  excitation  levels,  to  B  *  «,  corresponding  to  an  infinite  exci¬ 
tation  frequency.  In  the  case  B  *  0,  the  ionization  frequency  will  be  just  the 
inverse  of  the  time  required  for  the  electron  to  gain  an  amount  of  energy 
equal  to  the  ionization  energy. 


<V->  m  — 

1  €t 


hence. 
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(49) 


Ac/k 


As  a  first  approximation  it  will  be  assumed  that  B  is  small  and  that  the  growth 
rate  is  adequately  described  by 


v  A(/k 
x 


(50) 


The  electrons  in  the  laser  field  will  gain  energy  at  some  distribution 
varying  about  the  mean  rate  of  energy  gain.  A  general  distribution  function 
of  the  form 

j  (k>  dk  =  ^n- 


)' 


1 

r(t>) 


k*-  V  *k/kc 


dk 


(51) 


where  t)  is  a  numerical  factor  and  k^  is  the  mean  rate  of  energy  gain  will  be 
considered.  It  should  be  noted  that  equation  (51)  satisfies  the  conditions 


/* 

Jo 


j  (k)  dk  =  1 


(k)  dk  -  k 


(52) 


3 

If  a  value  t}  =  j  is  used  equation  (51)  has  a  Maxwellian  shape  about  the  mean 
energy  gain. 

The  mean  growth  rate  is 
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4 


Thua,  one  can  utilize  equation  (57)  to  study  the  average  electron  growth  rate 
for  a  very  generalized  distribution  in  the  rate  of  energy  gain. 

For  the  present  study  a  Maxwellian  type  of  distribution  function  will 
be  assumed,  thus, 


where 

z  kj(*)  =  *^5/2^  =  J  ®  (1  +  3z  *+  3z  2) 

and 


(59) 


(60) 


(61) 


It  should  be  noted  that  when  z  is  large  enough  so  that  3z”  1  and  3z"  2 
can  be  neglected  in  comparison  to  unity  then  equation  (59)  has  the  particularly 
simple  form 

The  only  parameter  in  equations  (59)  and  (62)  which  involves  the  laser  proper¬ 
ties  is  z,  hence,  the  relation  of  the  threshold  field  to  the  state  of  the  propaga¬ 
ting  medium  can  be  seen  quite  clearly. 

The  rate  at  which  an  electron  gains  energy  in  the  field  of  a  laser  beam 
is  given  by  the  relation 
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vibrational  levels  is  assumed  to  have  a  constant  cross  section  of 
-16  2 

3  x  10  cm  ,  over  the  energy  interval  1.  9  to  2.8  ev. 

For  the  singlet  transitions)  the  cross  section  is  given  by  equation  (14). 
Integrating  over  the  energy  Interval) 


cxA« 


■  •-  i'm- 

x  ' 


(66) 


(singlet  transition) 

The  integral  for  the  triplet  transitions)  using  the  cross  section  given 
by  equation  (19)> 


ox  Ac  =0.461 


</ 

X  €. 


Ci 


(c  -  C )  e 


* 


<U 


=  0.  17  y-  ^1  -  2.  72  —  e 


(67) 


(triplet  transitions) 


The  contributions  from  each  of  the  levels  is  given  in  Table  II. 
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Species ' 

State 

a  AG 

X 

(cm-ev) 

3  + 

A  I 

u 

6. 38  x  10"16 

N2 

C\ 

5.89  x  10‘17 

1 

a  n 

g 

3.  12  x  10‘16 

VIB. 

2.  70  x  10"16 

°2 

bV 

u 

3.42  x  10"16 

Table  U 

Contributions  From  Each  Level 


Assuming  a  composition  of  78%  N2  and  21%  02>  the  total  integral  of  the  exci¬ 
tation  cross  section  over  the  energy  range  from  threshold  to  ionization  is, 

—  - 15  2 

axAc  =  1*07  x  10  cm  -  ev  (68) 

The  contribution  from  the  N2  <c\  )  and  O^B^E^)  (Schumann-Runge  band) 
bands  is  negligible  in  comparison  to  the  other  states. 

The  calculated  results  for  the  complete  breakdown  of  air,  assuming 

3 

an  initial  electron  density  of  one  electron  per  cm  ,  are  shown  in  figure  9. 

The  results  of  calculations  by  Canavan  et  al.  (Ref.  16)  are  shown  for 

purposes  of  comparison.  The  agreement  is  relatively  good  but  is  somewhat 

9  2 

higher  than  the  quoted  experimental  value  of  p  ~  2  x  10  watts/cm  . 
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Nevertheless,  it  is  clear  that  the  critical  parameter  in  the  determin¬ 
ation  of  breakdown  for  a  given  laser  frequency  and  intensity  is  the  ratio 


vxA< 


(69) 


An  increase  in  vc  corresponds  to  an  increase  in  the  rate  at  which  energy  is 
absorbed  by  the  electrons.  In  the  calculations  carried  out  thus  far  it  has 
been  assumed  that  vc  was  exclusively  the  interaction  of  electrons  with  neu¬ 
tral  particles.  If  a  significant  degree  of  ionization  existed  at  the  time  the 
energy  was  being  absorbed,  this  would  be  written 


where  v  is  the  total  collision  frequency  between  electrons  and  neutrals  and 
cn 

v  is  the  total  collision  frequency  between  electrons  and  ions.  Because  the 

Cl  ^ 

-ross  section  for  collisions  between  electrons  and  ions,  interacting  through 
the  Coulomb  potential,  is  so  much  larger  than  the  electron-neutral  cross 
section,  as  soon  as  there  is  a  significant  degree  of  ionization  the  rate  at 
which  energy  is  gained  increases  rapidly  and  the  air  will  break  down  at  pro¬ 
portionately  lower  power  levels. 

The  implication  of  this  is  to  look  for  some  way  in  which  the  degree 
of  ionization  m.  the  air  is  increased  at  relatively  low  power  levels.  This 
usually  implies  the  presence  of  impurity  particles,  thus  one  must  look  at 
the  properties  of  dust  particles  in  air. 


Consider  a  particle  in  the  atmosphere  which  is  absorbing  radiation. 
The  change  in  the  temperature  at  the  surface  of  the  particle  will  be 


at 


SSL  At 

V 


(71) 


32 


where  Kie  the  absorption  coefficient,  <p  is  the  radiation  flux,  c  is  the 

v 

specific  heat  capacity  of  the  absorbing  particle,  and  p  is  the  particle  density. 
As  a  typical  example,  one  can  consider  quartz  (Si02)  which  has  an  absorp¬ 
tion  coefficient  of  K  ~  104cm-1  at  \  *  10.  64  ,  density  p«2.7  gm/cm3.  and 
specific  heat  capacity  cy  ~  1.  3  x  107  ergs/gm  *K.  The  vaporization  temper¬ 
ature  for  SiO^  is  T  -  2.  5  x  10  *K,  thus,  in  order  to  raise  the  temperature 

of  the  surface  of  the  particle  to  the  vaporization  level  in  a  time  t  ~10~®  sec. 

7  > 

requires  a  flux  <p  x  10  watts/cm  .  The  equilibrium  electron  density  for 
Si02,  which  exists  as  SiO  in  the  vapor  state,  can  be  obtained  from  the  Saha 
equation, 


where  h  is  the  Planck  constant  and  g  is  the  degeneracy.  Solving  for  N  , 

© 


where  S(T)  is  just  the  right  hand  side  of  equation  (72).  Taking  a  particle 
density  of  x  10  cm  3  and  an  ionization  energy  of  (  10.  8  ev,  the 

equilibrium  electron  density  at  vaporization  temperature  will  be  N  ~10llcm 
While  this  is  a  significant  increase  in  electron  density,  it  is  probably  not 
enough  to  lead  to  a  major  change  in  the  breakdown  power  levels.  Thus,  one 
must  do  a  somewhat  more  careful  study  of  the  heating  and  subsequent  vapor¬ 
ization  of  dust  particles  with  particular  attention  to  various  mechanisms  for 
electron  enhancement. 
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SECTION  HI 

THE  HYDRODYNAMIC  EQUATIONS 

*  ‘  n  1  l  ,  ,  ;  .  . 

- 

Two  separate  models,  one  utilising  a  single  temperature  for  both 
ions  and  electrons  and  the  other  model  treating  the  temperatures  separately, 
have  been  developed  for  the  description  of  dirty  air  breakdown.  For  this 
report,  however,  only  the  more  realistic  two -temperature  model  has  been 
used  for  calculations.  In  both  models  a  symmetry  is  assumed  so  that  the 
problem  can  be  treated  in  one  dimension.  The  Lagrangean  form  of  the  one¬ 
dimensional  hydrodynamic  equations  is  used.  The  first  three  equations  are 

both  models  (Ref.  17), 

=  U  (74) 


essentially  the  same  in 


dR 

dt 


au 

dt 


R“'1  a!i<p+Q> 


(75) 


V=  1  =  1 

p  a  d  M 


(76) 


expressing  the  definition  of  velocity  U,  the  conservation  of  momentum  and  the 
conservation  of  mass.  The  artificial  viscosity  Q  is  necessary  for  numerical 
stability  of  the  difference  equations  discussed  in  Appendix  I.  The  geometry 
of  the  single  dimension  is  plane,  cylindrical  or  spherical  for  a  s  1,  2  or  3, 
respectively,  and  M  is  the  Lagrangean  mass  coordinate. 

The  one-temperature  model  is  based  on  the  assumption  that  the  free 
electrons  are  in  thermal  equilibrium  with  the  neutral  molecules'  and  ions  of 
the  shock  heated  air  and  blow  off  vapor.  The  electron  density  is  determined 
by  the  Saha  equation.  The  energy  equation  of  hydrodynamics  is 
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2E 

at 


(77) 


_  p  il+  Tr  I  5  f  Da- 1  ST  1 

p  ir+  *K+  ns  I*r  ir  I 

where  x  *■  the  thermal  conductivity,  p  it  the  laser  flux  and  K  is  the  effec-  r' 
tive  mass  absorption  coefficient  for  the  material  and  free  electrons.  The 
hydrodynamic  description  is  completed  by  air  equation  of  state  giving  the 
pressure  P  as  a  function  of  the  density  pand  specific  internal  energy  E. 

.  In  the  two-temperature  model,  it  is  assumed  that  the  electrons  are  in 

a  Boltzmann  distribution  at  temperature  T  .  The  electrons  interact  with  the 

£ 

neutral  molecules  and  ions  which  are  in  a  second  Boltzmann  distribution 
with  temperature  T^.  The  total  pressure  in  the  momentum  equation  is  the  ' 
sum  of  the  partial  pressures  of  the  electrons  and  ions. 

P  s  Pe  +  Pi  •  (78) 


For  notation,  the  neutral  molecules  and  ions  are  referred  to  collectively  as 
"ions.  " 

The  energy  equation  qow  has  separate  components  for  the  ions  and 
electrons. 


dE 
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where  IJ+P^sP+Q,  v  represents  a  collision  frequency,  c  is  an  energy, 
and  E  is  a  specific  energy.  The  subscripts  c,  x,  a  and  i  represent  elastic 
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collisions,  excitation  collisions,  attachment,  and  ionisation,  respectively. 
The  change  in  the  electron  density  is  given  by 


r* 


+  V.  N 
i  e 


v  N 

a  e 


(81) 


where  is  the  electron  diffusion  coefficient.  Functional  forms  for  these 
parameters  are  discussed  in  the  next  section. 


1.  COLLISION  FREQUENCIES 


The  energy  losses  within  the  electron  and  ion  gases  are  expressed  in 
terms  of  the  various  collision  frequencies  averaged  over  some  suitable  dis¬ 
tribution  of  energies.  For  the  calculations  to  be  carried  out  in  this  report, 
it  will  be  assumed  that  the  components  of  the  gas  have  a  Maxwellian  velocity 
distribution  about  some  temperature  representative  of  the  component  under 
consideration,  thus 


(82) 


where  N  is  the  density  of  target  particles,  m  is  the  mass  of  the  colliding 
particle  (in  some  cases  the  reduced  mass),  and  a  is  the  cross  section  for 
the  interaction  under  consideration.  The  integration  is  carried  out  over  the 
appropriate  energy  range. 


a.  Electronic  Excitation 

Utilizing  equation  (14)  of  Section  II  for  the  cross  section  for 
singlet  transitions,. 


e-c/*Tdf 


(83) 
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where  m  ie  the  electron  mast,  e  is  the  electron  charge,  and  T  is  the 
electron  temperature.  The  integration  is  straightforward. 


(singlet  transition) 


where  E^z)  is  an  exponential  integral  of  the  first  kind.  In  the  special  case 


C  >  >  kT  , 
x  e  . 
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(singlet  transition) 

A  similar  calculation  for  the  triplet  transitions  gives 


(86) 


0.  340  e4N 
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(triplet  transition) 


b.  Vibrational  Excitation 

The  cross  section  for  the  resonant  excitation  of  N.  vibrational 

€m 

levels  has  been  assumed  to  be  constant  over  the  energy  range  of  interest, 
hence, 


VkT, 


(87) 
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where  oy  m  3  x  10  cm  . ,  (,  <*»  1.8  ev  and  *2.  8  ev. 


c.  Ionization  Frequencies 

For  N2,  using  equation  (25)  for  the  cross  section 


-15  2 

1.64  x  10  V 


4  -<i/kT, 


and  for  O^,  using  equation  (26), 


1.74  x  10"15f.2 


Electron-Ion  Elastic  Collision  Freauenc 


The  elastic  collision  frequency  for  two  particles  interacting 
through  a  Coulomb  potential  can  be  evaluated  from  classical  considerations 
(Ref.  18).  The  result,  after  averaging  over  the  relative  velocities  is. 


(90) 
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where  m  is  the  reduced  mass  and 
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The  length,  rD  ,  is  the  well  known  Debye  length  and  N.  ,  N  are  the  ion  and 
electron  number  densities,  respectively.  For  the  case  of  electron-ion 
collisions,  the  reduced  mass  in  equation  (90)  is  replaced  by  the  electron 
mass. 

2.  AMBIPOLAR  DIFFUSION  COEFFICIENT 

When  the  Debye  length  in  a  plasma  is  small  in  comparison  to  the 
characteristic  diffusion  length,  the  free  electron  diffusion  coefficient  must 
be  replaced  by  the  ambipolar  diffusion  coefficient.  The  Debye  length  ii 
defined  by 


For  typical  conditions  of  interest,  T  3  x  103  *K,  N  >  1015cm'3.  the 

-5  e 

Debye  length  is  r_  10  cm.  The  free  electron  diffusion  coefficient  under 
U  2  2 

the  same  conditions  will  be  D  «  10  cm  / sec.  The  characteristic  relaxa- 

e  -8 

tion  time  of  interest  will  be  <r~  10  sec. ,  hence, 
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cm. 


(94) 


L_  **  \/  D  t'  ~  1 0 
D  •  *  e 

Thus,  cne  must  consider  the  space  charge  effects  on  the  diffusion  of  the 
electrons. 

To  do  so  one  notes  that  the  effect  of  the  space  charge  will  be  to  cause 
the  electron  to  pull  an  ion  along  with  it.  This  can  be  viewed  as  increasing 
the  mass  of  the  electron  in  the  electron  diffusion  coefficient  by  the  mass  of 
the  ion.  The  ion  collision  cross  section  with  heavy  particles  will  be  about 
twice  that  of  the  electron  cross  section,  hence, 


(95) 


This  is  correct  only  when  the  ion  and  electron  temperatures  are  equal. 
However,  for  purposes  of  the  present  report  the  temperatures  will  be  close 
enough  over  times  of  interest  so  that  equation  (95)  can  be  used. 


3.  FREE -FREE  ABSORPTION  COEFFICIENT 

The  interaction  between  a  current  and  an  electric  field  is  described 
by  the  Maxwell's  equations, 


Reducing  equations  (96)  and  (97)  to  a  single  wave  equation 


(96) 

(97) 
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(98) 


+ 


4  ff  &£ 
2  dt 


where  "j  is  the  current  density, 

?=-Nev  (99) 

where  N  is  the  electron  density  and  v  is  the  electron  velocity. 

The  motion  of  an  electron  is  an  electric  field  is  described  by  the  equa¬ 
tion  of  motion 


dv  .  ei£ 

TT  +  V  v  - - 

dt  c  m 


(100) 


where  v.  is  the  collision  frequency.  Assuming  that  the  vector  quantities  all 
have  the  time  behavior  e"1  wt 
equations  (98)  and  (100) 


where  u  is  the  oscillation  frequency,  reduces 


v  =  - 


ei£  (v  +  i u>) 

c _ 

m(w2+  v2) 


(101) 


c2  c2  ‘ 


Substitution  of  equation  (101)  into  equation  (102) 

72^  +  k2£  =  0 


(102) 


where 
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4rre  N(v  +  iw) 

k  =  £  1  +  i - — 

c  ,  2,  2. 

mo)(o)  +  v  ) 


(103) 


Assuming  that  the  second  term  in  the  square  brackets  is  small  in  compari¬ 
son  to  unity,  one  can  use  the  usual  approximation. 


2  ffe  N(v  +  i w) 
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c  .11 
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(104) 


The  absorption  coefficient  will  be  just  the  imaginary  portion  of  equation 
(104) 


K„  = 
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This  is  the  usual  form  for  the  free-free  absorption  coefficient.  The  mass 

■A 

absorption  coefficient  of  equation  (105)  is  just 


-  Kff 

K  =  — =- 
e  p 


(106) 
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SECTION  IV 


PARTICLE  DESCRIPTION  AND  INTERACTIONS 

The  size  of  particles  in  dry  atmospheric  aersols  range  from  about 
0.005^  to  20|i.  Smaller  particles  are  quickly  attached  to  larger  particles 
and  are  found  only  in  a  region  of  constant  production.  Particles  larger 
than  about  20(i  are  lost  by  sedimentation  and  are  found  only  in  the  vicinity 
of  a  source.  The  density  of  particles  versus  radius  shows  a  maximum 
near  10^  (Ref.  19). 

Air  breakdown  experiments  have  been  performed  on  particles  in  the 
range  40-l00|j  (Ref.  2).  In  order  to  compare  our  results  with  these  ex¬ 
periments.  computer  runs  were  made  for  particles  ranging  in  diameter 

•.  ,  * 

from  12|im  to  96|im. 

1.  ENERGY  ABSORPTION  FROM  THE  LASER  BEAM 

Energy  is  absorbed  from  the  laser  beam  by  two  processes:  mass 
absorption  by  the  particle  material,  and  inverse  bremsstrahlung  absorption 
r  t  fey  ^ee  electrons.  As  the  particles  are  heated  and  ionized,  the  mechanism 
for  absorption  of  energy  from  the  laser  beam  shifts  from  the  first  to  the 
second  kind. 

a.  Mass  Absorption 

In  general,  it  is  necessary  to  know  the  mass  absorption  coeffi¬ 
cients  of  the  various  types  of  dust  particles  which  may  be  encountered  in 

the  atmosphere.  Specifically  under  consideration  are  the  materials  A1  O  , 

2  3 

SiO^.  C.  and  NaCl  as  representative  of  the  common  types  of  particles 
found  in  the  atmosphere.  The  absorption  coefficient  is  needed  as  a  function 
of  wavelength  (at  least  for  specific  laser  wavelengths  of  interest),  as  a 
function  of  particle  material  temperature,  and  as  a  function  of  material 
density  as  the  particle  melts  and  vaporizes. 
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Data  of  the  kind  needed  are  scarce.  Data  are  available  for  normal 
density  solids  in  the  wavelength  range  of  particular  interest  in  optical 
(Ref.  20)  and  infrared  (Ref.  21 )  applications  at  room  temperatures.  Good 
data  for  solids  at  elevated  temperatures  are  known  only  for  sapphire  (Ref. 

22)  in  the  0.  5  to  6^m  wavelength  range.  An  analysis  of  the  absorption  co¬ 
efficient  of  quartz  vapor  (SiO)  in  air  has  been  give  i  for  10.  6um  wavelength 
as  a  function  of  temperature  and  density  (Ref.  23).  The  reflection  coeffi¬ 
cient  for  neodymium  laser  light  on  a  carbon  target  and  targets  of  otner 
materials  has  been  reported  as  a  function  of  laser  pulse  power  (Ref.  24) 
level.  Common  salt  is  quite  transparent  at  most  wavelengths  of  interest. 

In  all,  the  available  empirical  data  are  far  from  complete,  and  further  ex¬ 
perimental  measurements  in  this  area  would  be  highly  desirable.  The  known 
information  is  summarized  below. 

(1)  sio2 

Extensive  data  exist  (Refs.  20  and  25)  on  the  values  of  the 
absorption  coefficient  of  Si02  in  its  various  forms  at  room  temperature  for 
wavelengths  in  the  range  from  0.2  to  4.  8|jm.  The  temperature  dependence 
is  not  known  to  have  been  measured,  however,  as  it  has  for  sapphire.  When 
Si°2  *■  heated  in  air,  it  forms  SiO  vapor  as  a  result  of  oxidation.  The  ex¬ 
tinction  coefficient  of  SiO  has  been  measured  (Ref.  26)  as  a  function  of  wave* 
length  from  0.  24  to  I4|im.  A  plot  of  the  absorption  coefficient  from  these 
data  is  shown  in  figure  11  for  SiO.  Also  shown  in  the  figure  is  the  room 
temperature  absorption  coefficient  for  silica,  and  one  point  measured  at 
AFWL  for  fused  silica  at  2.8^m.  A  value  of  the  absorption  coefficient  be¬ 
tween  2000*  and  about  4000*K  for  normal  density  (p  =  2.  13  gm/cc)  SiO 
derived  from  the  work  of  Chang,  Drummond,  and  Hall  is  also  shown  for 
comparison.  It  falls  below  that  for  the  room  temperature  material  despite 
being  for  elevated  temperatures.  Both  values  were  used  in  the  computer 
runs  made. 
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The  density  dependence  of  the  absorption  coefficient  employed  was 
of  the  form  (p7 pQ)  ,  where  n  was  chosen  to  equal  two  for  SiO  in  conform¬ 
ance  with  the  value  given  for  diatomic  molecules  by  Chang,  Drummond, 
and  Hall. 

(2)  A1203 

Measured  values  of  the  absorption  coefficient  for  sapphire 
as  a  function  of  wav&length  and  temperature  are  shown  in  figure  10.  Most 
of  these  data  comes  from  the  work  of  Gryvnak  and  Burch  (Ref.  22).  A  few 
points  reported  in  1936  are  also  shown  but  are  of  questionable  accuracy  be¬ 
cause  of  suspected  impurities  (Ref.  27).  Two  additional  points  are  shown. 
One  was  measured  at  AFWL  (Ref.  28 );  the  other  reported  as  the  limit  of 
the  transmission  band  in  reference  21.  These  individual  points  were  all 
for  room  temperature  material.  Gryvnak  and  Burch  (Ref.  22)  also  report 
a  qualitative  estimate  that  the  absorption  coefficient  increases  discontinu¬ 
ous^  by  a  factor  of  the  order  of  30  when  sapphire  melts  (at  2045*C).  As 
a  first  approximation,  it  may  be  assumed  tiat  this  value  remains  about 
constant  up  to  the  vaporization  temperature  (2980 *C). 

Extending  the  known  data  to  cover  the  complete  range  of  interest  is 
largely  speculative.  The  trend  of  the  known  curves  in  figure  10  have  been 
extrapolated  af[  shown  to  the  C02  laser  wavelength  of  interest.  A  value  of 
about  l.Oum  is  estimated  on  this  basis  for  the  absorption  coefficient  of 
sapphire  at  10.  6  4m  wavelengths  at  room  temperature.  The  temperature 
dependence  at  5  4m  wavelength  is  roughly  linear  on  a  log -log  plot,  so  an 
empirical  formula  for  the  below r melt  temperature  dependence  factor  is 
suggested  of  the  approximate  form. 


f(T) 


T  -  300 
790 


for  T  <  2318*K 


(107) 


where  T  is  the  temperature  in  degrees  Kelvin  and  f(T)  is  a  multiplying 
factor  for  the  room  temperature  absorption  coefficient.  Above  ?.318#K  a 
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Figure  10.  '  Absorption  Coefficient  for  A^O^  as  a 
Function  of  Wavelength 
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constant  factor  of  the  order  of  400  can  be  assumed  to  multiply  the  room 
temperature  value  of  the  absorption  coefficient. 


In  conformity  with  the  findings  of  Chang,  Drummond  and  Hall  (Ref. 
23),  the  dependence  of  the  absorption  coefficient  of  sapphire  vapor  can  be 
further  taken  to  depend  on  density  as  some  power  of  the  density. 


.n 


•  Chemical  reactions  will  take  place  when  Al^O^  is  vaporized  in  air 

which  alter  its  composition.  A  detailed  theoretical  study  of  this  problem 
was  considered  beyond  the  scope  of  this  effort.  For  purposes  of  estima¬ 
tion,  an  absorption  coefficient  was  essumed  which  was  independent  of  tern* 

perature  but  depended  on  the  density  as  (p/p J2  in  conformity  with  the  one 

c 

employed  for  silica  in  the  computer  runs  made. 


(3)  NaCl 

The  room  temperature  absorption  coefficient  for  NaCl 

is  shown  in  figure  12  at  wavelengths  above  10  um.  At  the  CO  laser  wave- 

£ 

length  of  10.  6  urn,  a  small  particle  of  NaCl  is  effectively  transparent,  and 
it  remains  so  down  to  about  0.  2  ^im.  Consequently,  at  the  laser  powers  of 
interest  in  this  study,  the  heating  would  be  very  slow  relative  to  that  of  the 
other  materials  investigated  at  10.  6  (am.  No  data  are  known  for  NaCl  at 
elevated  temperatures. 

(4)  Carbon 

Carbon  takes  on  many  forms  which  vary  greatly  in  their 
optical  properties.  Not  only  is  diamond  transparent  from  0.  3  to  80  ym, 
but  very  pure  crystals  of  graphite  have  also  been  found  which  are  trans¬ 
parent  (Ref.  28)  in  the  visible  region.  On  the  other  hand,  usual  forms  of 
carbon  are  opaque.  One  reason  given  for  this  is  the  density  of  dislocations 
in  graphite  crystals  determines  the  ease  with  which  light  can  pass  through 
them.  Another  reason  given  is  that  most  carbon  is  very  porous,  and  the 
light  entering  it  is  multiply  scattered  before  it  can  escape  »*>d  is  thus  ab¬ 
sorbed.  Soot,  which  may  have  a  density  two  order  of  magnitude  below 
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solid  graphite  crystals,  is  particularly  porous,  and  as  a  consequence  is 
used  as  an  antireflecting  coating  in  the  form  of  lampblack.  A  study  has 
been  made  of  the  transmittance  of  thin  layers  of  soot  and  other  carbon 
blacks  (Ref.  29).  From  these  data,  an  effective  mass  absorption  coeffi¬ 
cient  for  soot  like  carbon  may  u.c  calculated,  and  the  results  are  shown  in 
figure  *13. 

The  effect  of  temperature  on  the  absorptivity  of  carbon  is  not  known 
to  have  been  measured  directly.  However,  that  it  increases  may  be  infer¬ 
red  from  the  observed  decrease  in  the  reflection  coefficient  for  light  from 
intense  neodymium  laser  pulses  as  the  fluence  was  increased  (Ref.  24). 
Because  of  the  uncertainty  in  the  form  of  carbon  particles  which  might  be 
encountered  in  the  atmosphere,  and  because  no  experiments  were  known 
which  used  carbon  particles  to  reduce  breakdown,  further  work  on  the  ab¬ 
sorptivity  of  carbon  was  discontinued. 

For  all  particle  materials,  the  vapor  phase  tends  to  exhibit  a  reduced 
absorption  as  the  band  structure  of  the  solid  is  lost.  In  the  limit  of  zero 
density,  the  mass  absorption  coefficient  goes  to  zero,  except  in  the  unlikely 
case  that  one  of  the  discrete  lines  of  the  molecular  absorption  spectra  is 
exactly  at  the  laser  frequency.  At  higher  density,  collision  broadening  of 
the  molecular  levels  can  cause  an  increase  in  absorption.  The  mass  absorp¬ 
tion  coefficient  is  given  by 


A  =  -  (108) 

P 

where  a  is  the  linear  absorption  coefficient.  The  density  dependence  of  the 
absorption  is  expressed  in  the  form 


a(p)  *  a 


0 


(109) 
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Wavelength  -  (microns) 


Mass  Absorption  Coefficient  for  Carbon 
(soot)  as  a  Function  of  Wavelength 
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where  a^  is  the  absorptivity  of  the  solid,  is  the  density  of  the  solid,  and 
n  is  to  be  determined.  For  X-ray  deposition,  the  energy  is  deposited  in 
;he  atomic  electrons,  n  =  1  and  the  mass  absorption  is  independent  of  den¬ 
sity.  For  laser  deposition  in  the  vapor  from  a  quartz  particle  n  -  2  (Ref. 
23).  For  this  diatomic  molecule,  one  density  factor  is  from  the  density  of 
absorbing  molecules  and  one  from  the  collision  frequency  which  causes  line 
broadening.  Numerical  studies  have  shown  that  the  hydrodynamic  results 
are  insensitive  to  the  choice  of  this  parameter  for  n  a  1.  5  (Ref.  30).  Con¬ 
sequently,  we  have  used  the  value  n  =  2  for  all  materials  considered. 

b.  Absorption  by  Electrons 

Independent  of  the  mass  absorption  of  the  material  particle  and 
vapor,  free  electrons  can  absorb  energy  directly  from  the  laser  beam  by 
inverse  bremsstrahlung.  Absorption  by  this  mechanism  is  discussed  in 
Section  III- 3. 

2.  EQUATION  OF  STATE 

To  complete  the  hydrodynamic  description  of  the  system,  the  equa¬ 
tion  of  state  must  be  specified.  The  material  particle  is  described  by  an 
equation  of  state  of  the  Mie-Grilneisen  form  (Ref.  31) 

P*P0+rp(E-E0)  (110) 

where  r  !•  the  Grttneisen  coefficient.  Pq(p)  and  Eg(p)  give  the  density 
dependence  of  the  pressure  and  energy  density  at  a  reference  temperature. 
In  the  comprassional  region, 

u  =  -2-  -  1  >0  (111) 


the  pressure  is  given  by 
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p  =  PH(u>  [  1 


(112) 


where 


Pjj(|i)  =  Cji  +  Du2  +  S^3 


(113) 


and  the  Hugoniot  coefficients  are  determined  from  experimental  data  (Ref. 
32).  In  the  expansion  region, 


(1U) 


P  = 


♦(T-yj)* 


where  Eg  is  the  sublimation  energy, 


(115) 


Yi  =  r- 1 


(116) 


and  y  is  the  ratio  of  specific  heats  for  the  vapor.  The  factor 


F  * 


rE, 


SPc 


ensures  continuity  with  the  compress ional  equation  as  **  •*  1. 


(117) 


The  expansion  equation  must  also  allow  for  a  change  of  phase  as  the 
material  is  heated.  In  the  above  expression,  the  phase  transition  is  included 
only  in  an  average  or  smoothed  way  as  discussed  in  Ref.  32.  Near  normal 
density  and  energy  density,  the  expansion  equation  gives  a  pressure  appro¬ 
priate  to  the  solid.  At  low  densities  and  high-energy  density,  the  expansion 
pressure  approaches  that  of  an  ideal  gas.  For  the  problem  of  dirty  air 
breakdown,  the  region  near  the  phase  transition  is  critical  for  it  controls 


53 


the  rate  of  material  vaporization,  blowoff  and  aubsequent  shock  heating. 
There  is  experimental  evidence  (Ref.  33)  that  for  some  materials  the  blow- 
off  consists  of  macromolecular  complexes  rather  than  single  molecules. 
This  would  completely  change  the  above  model;  changing  the  effective  sub¬ 
limation  energy  of  the  solid  and  tr«»  y  of  the  vapor. 

For  the  air  surrounding  the  dust  particle,  too  separate  models  have 
been  used.  In  the  first  model,  the  air  is  assumed  to  be  in  local  thermo¬ 
dynamic  equilibrium,  including  the  various  molecular  processes.  The 
pressure,  equilibrium  electron  density  and  effective  specific  heat  ratio  for 
shocks  are  computed  as  functions  of  density  and  energy  density  using  an 
interpolation  scheme  on  tables  given  by  Gilmore  (Ref.  34). 

The  validity  of  the  equilibrium  assumption  is  questionable  since  the 
rate  constants  for  the  molecular  processes  are  often  long  compared  to  laser 
pulse  times  of  a  few  microseconds.  For  example,  at  standard  atmospheric 
density  and  3000 *K,  the  reaction 

N2  +  °2  *  2NO  (118) 

requires  about  10~4  second  to  reach  equilibrium  (Ref.  35).  This  particular 
reaction  is  important  for  the  equilibrium  electron  density  because  the  ion¬ 
ization  potential  of  NO  at  9-  25  ev  (Ref.  36)  is  the  lowest  of  the  molecular 
species  found  in  air. 

For  laser  pulses  that  are  short  compared  to  the  molecular  reaction 
times,  a  better  approximation  to  the  air  equation  of  state  is  obtained  by 
assuming  no  molecular  dissociation  or  reaction.  This  is  the  basis  of  the 
second  model  in  which  the  pressure  is  given  by  the  ideal  gas  equation.  The 
equilibrium  electron  density  is  computed  using  the  Saha  equation  with  a 

single  ionization  potential  of  12.  05  ev  for  O  . 

2 

Calculations  were  carried  out  using  both  models.  Over  the  times  of 
the  order  of  10  sec.  there  was  no  significant  difference  in  the  answers 
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obtained  from  the  two  models.  The  results  presented  in  later  sections  were 
all  obtained  from  the  second  model. 

The  equation  of  state  for  the  electron  gas  is  just  that  of  an  ideal 
monatomic  gas. 


Pe  =  (ye-i)pEe  (119) 

5 

where  y  =  —  and  E  is  the  electron  internal  energy  per  gram  of  material. 

©  <3  © 

Properties  of  the  materials  studied  art  listed  in  Table  Ill.  The  para¬ 
meter  values  employed  in  the  CELAB  computer  runs  are  given. 

3.  THERMAL  CONDUCTIVITY 

In  general,  the  thermal  conductivity  of  the  dust  particle  will  depend 
on  the  temperature.  We  have  assumed  the  functional  form 

X  (T)  =  Xq  T^  (120) 

for  the  temperature  dependence.  We  have  fit  limited  experimental  data  on 
quartz  and  sapphire  to  determine  values  of  and  for  the  two  materials. 

A  plot  of  the  fitted  conductivity  is  shown  in  figure  14. 

The  thermal  conductivity  of  air  shows  a  general  increase  from  room 
temperature  to  about  7000*K,  then  decreases  between  7000  and  10,  000*K 
as  shown  in  figure  15.  The  data  are  too  complicated  to  be  fit  with  the  func¬ 
tional  form  of  equation  (120)  with  a  single  value  of  the  parameters  over  the 
entire  temperature  range.  We  have  divided  the  data  into  four  separate 
regions  and  fit  each  region  with  the  form  of  equation  (120).  For  any  given 
temperature,  a  table  look-up  is  used  to  determine  the  values  of  the  coeffi¬ 
cients  and  the  conductivity  computed  from  the  above  form. 

The  thermal  conductivity  of  the  free  electron  gas  can  be  computed 
from  the  kinetic  theory  of  gases  (Ref.  37). 
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TABLE  IU 


MATERIAL  CONSTANTS 


Property 

Air 

Al2°3 

Density  (gm/cc) 

1.292  x  10“3 

2.  65 

3.97 

y 

1.  3802 

1.25 

1.25 

Atomic  Weight 

28.  98 

60.06 

102.0 

Thermal  Conductivity  Constant 
(ergs/cm  sec"K^+* ) 

31.4 

1-820 

1820 

Linear  Absorption  Coeff  (cm  ) 

0.78 

0. 

0.74 

1  X  io4 

0.74 

1  x  io4 

GrUneisen  Coefficient 

:  --- 

0.78 

0.78 

Hugoniot  Constants 

8.7  x  10H 

4.  65  x  1012 

2 

C  (dyne/cm  ) 

— 

2 

D  (dyne/cm  ) 

0 

0 

S  (dyne /cm2) 

— 

0 

0 

Sublimation  Energy  (erg/gm) 

— 

2.  35xlOU 

2.  35  X  10H 

Specific  Heat  Solid  (erg/gm  deg) 

— 

1.25X107 

1. 72  x  10? 

Vaporization  Temp  (*K) 

— 

2503 

3253 

Work  Function  (ev) 

— 

4.  58 

3.77 

Ionization  Energy  (ev) 

Slope  of  Ionization  Cross 

12.  05 

4.  5  x  10"2 

10.8 

8. 3  x IO-2 

7.7 

8i3x  10"2 

2 

Section  (ffa^  /ev) 

Vibration  Threshold  (ev) 

1.8 

1.8 

\ 

1.8 

Vibration  Bandwidth  (ev) 

1.4 

1.4 

1.4 

Vibration  Quanta  (ev) 

0.26 

0.26 

0.26 

Vibrational  Cross  Section 

,  n- 1 5  2. 

(1  x  10  cm  ) 

26.  60 

26.  60 

26.  60 

K(10  erg/cm* 


T(103,K) 


Figure  15.  Air  Conductivity,  K,  and  Diffuiivity, 
D,  at  One  Atmosphere  Pressure. 
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x0=fN,kD«  <u‘> 

where 

De  =  }  v  Xe  (122) 

is  the  electron  diffusion  coefficient,  v  is  the  average  electron  velocity  and 

\  is  the  mean  free  path.  In  terms  of  the  collision  frequency,  the  temper* 
c 

ature -dependent  diffusion  coefficient  is  given  by 


e  mv 

c 


Ambipolar  diffusion  is  discussed  in  Section  IQ-2. 


(123) 


4.  THERMIONIC  EMISSION 

The  thermionic  emission  (Ref.  38)  from  a  material  at  high  temper¬ 
atures  can.  be  evaluated  by  assuming  that  the  material  emitted  outside  of  the 
surface  is  in  equilibrium  with  the  material  inside.  Using  the  Fermi  distribu¬ 
tion  for  the  electrons, 


£(<  +  €k)  S  (\Y+<If)/kT  (124) 

1  +  e  K 

where  c  1*  the  energy  required  to  remove  an  electron  from  some  reference 
energy  level  to  a  state  of  rest  outside  the  metal,  <  is  the  kinetic  energy  of 
an  electron  outside  of  the  metal,  and  W  is  the  work  function.  Multiplying 
equation  (124)  by  N  v,  assuming  that  (W+fv-)  >  >  kT,  and  integrating  over  the 
density  of  states  in  phase  space  gives 
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J  =  — !py  (kT)2  c'W/kT 
2«  h 


,  ia20  2  -W/kT  -2  -1 

=  7.72  x  10  T  c  cm  -  sec 


(125) 


This  is  juat  the  well  known  Richardson  equation  for  thermionic  emission. 

On  the  basis  of  the  derivation  given  above,  one  expects  that  the  flux 
of  electrons  for  different  materials  would  depend  on  the  density  of  electrons 
in  the  conduction  band.  Thus,  to  some  degree,  the  difference  in  the  conduc¬ 
tivity  between  materials  should  indicate  the  difference  to  be  expected  in  the 

7 

thermionic  emission.  The  conductivity  of  most  metals  is  about  10  mhos/ 

•  IT 

meter  while  the  conductivity  for  quartz  is  ~  10  mhos/meter.  Thus,  one 
would  expect  the  thermionic  emission  for  quartz  to  be  down  several  orders 
of  magnitude  from  that  given  by  the  Richardson  equation. 

However,  experimental  evidence  (Ref.  39)  indicates  that  there  is  a 
substantial  flux  of  electrons  from  quartz  and  borosilicate  targets.  This  flux 
has  been  interpreted  as  being  due  to  heating  of  surface  contaminants.  This 
effect  would  show  up  in  the  Richardson  equation  as  a  change  in  the  constant 
factor.  Since  it  is  not  clear,  without  further  study,  what  value  the  constant 
factor  should  have  for  SiO^  and  Al^Oy  the  value  given  in  equation  (125)  has 
been  used  for  the  preliminary  calculations  given  in  this  report.  This  provides 
an  upper  limit  on  the  flux  of  electrons  due  to  thermionic  emission. 

5.  MACROMOLECULAR  COMPLEXES 

A  further  source  of  electrons  may  result  from  the  interaction  of  a 
laser  beam  with  naturally  occuring  molecular  clusters.  Such  clusters  gener¬ 
ally  have  a  substantially  lower  ionization  potential.  Of  primary  importance 
would  bo  water  clusters  in  the  atmosphere.  The  evaluation  of  the  density, 
size,  and  ionization  potential  of  such  clusters  is  a  difficult  problem  for  which 
there  exists  little  Information.  The  problem  is  neglected  in  the  present 


report. 
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SECTION  V 


RESULTS 

Calculations  have  been  carriod  out  to  determine  the  effect  of  impur¬ 
ities  on  the  breakdown  of  air  by  a  CO^  laser  beam.  The  particular  impur¬ 
ities  considered  were  sapphire  (Al^O^)  and  quartz  (Si02>,  although  the 
CELAB  hydrodynamic  code  developed  for  this  study  is  quite  general  and 
can  be  utilized  for  the  study  of  other  materials. 

For  numerical  considerations  in  the  present  report,  the  only  particle 
sizes  considered  were  12,  24,  48  and  96  ,  although  the  code  is  capable  of 

handling  particles  of  any  size.  The  code  divides  the  particles  into  a  number 
of  zones,  each  zone  having  a  width  of  approximately  10  y  .  In  air,  immed¬ 
iately  outside  the  particle,  the  zone  width  is  increased.  The  hydrodynamic 
time  steps  are  of  the  order  of  10  ^  seconds.  The  temperatures,  both 
ionic  and  electronic,  are  considered  throughout  the  particle  and  into  the  air. 

Figure  16  shows  the  calculated  electron  density  in  the  first  air  zone 

as  a  function  of  the  time  for  four  different  sizes  of  A1  O  in  the  field  of  a 

8  ^  *  2 

10.  6  4  laser  beam  having  an  intensity  of  3  x  10  watts/cm  .  The  peak 

electron  density  is  for  the  largest  particle  (96  y),  however,  the  difference 

in  the  electron  density  for  the  various  sizes  is  net  very  significant.  The 

electron  density  in  the  first  zone  reaches  the  steady  state  value  of  about 
16  3 

8  x  10  electrons/cm  .  This  represents  the  degree  of  ionization  of  about 
0.3%. 

A  similar  pair  of  curves  for  SiO,  and  sizes  24  and  48  (j  is  shown  in 

^  8  2 
figure  17.  Again,  the  laser  is  10.  6  n  with  an  intensity  of  3  x  10  watts/cm  . 

A  steady-state  value  of  the  electron  density  in  the  first  zone  is  apparently 

_7 

obtained  at  about  1.5  x  10  seconds. 

The  curves  for  SiO.,  figure  17,  show  a  pronounced  maximum  and 

L  _o 

minimum  occurring  at  approximately  3.  5  and  6  x  10  seconds.  Again,  the 
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Figure  17.  Electron  Denaity  va  Time  for  Different  Size  SiO,  Particle* 
in  the  Beam  CC>2  Laaer  at  an  Intenaity  of  3  x  10®  watta/cm^. 


difference  in  the  electron  density  for  the  two  sizes  is  not  substantial.  How¬ 
ever,  the  difference  in  the  peak  electron  density  for  SiO^  as  compared  to 
A1  O  is  about  4  orders  of  magnitude. 

The  maximum  and  minimum  of  the  curves  for  SiO.  indicate  a  diffusion 

L  -8 

type  of  motion  with  an  apparent  characteristic  time  of  about  3  x  10  seconds. 
The  curves  indicate  that  a  peak  density  builds  up  in  the  first  zone  at  some 
point,  shielding  the  remaining  portion  of  the  particle  from  further  heating  by 
the  laser  beam  aft  the  electrons  in  the  first  zone  absorb  most  of  the  energy. 
The  electrons  diffuse  out  of  the  zone  with  a  relatively  large  velocity.  When 
the  electron  loss  from  the  zone  approaches  some  limiting  value,  the  particle 
again  begins  to  absorb  most  of  the  laser  energy,  heating  up  and  emitting  more 
electrons  into  the  first  zone.  Since  the  diffusion  has  slowed  down,  the  second 
peak  is  more  persistent,  leading  to  a  quasi-steady  state. 

If  this  interpretation  is  correct,  it  would  indicate  a  diffusion  coeffi- 

2  2 

cient  of  approximately  10  cm  /seconds.  This  value  is  more  characteristic 
of  the  free  electron  diffusion  coefficient  than  the  ambipolar  diffusion  coeffi¬ 
cient  which  would  be  about  2  to  3  orders  of  magnitude  smaller.  Thus,  it 
appears  that  the  electrons  have  been  allowed  to  diffuse  out  of  the  zone  much 
too  rapidly. 

Figure  18  shows  the  electron  density  as  a  function  of  time  for  a  24  u 

8,2 

and  48  quartz  particle  at  a  CO  laser  intensity  of  10  watts/cm  .  These 

2  7-1 

curves  have  a  growth  rate  on  the  order  of  10  sec 

Figure  19  shows  an  electron  density  of  an  A1  O  particle  in  the  field 

7  3  7  2 

of  a  COz  laser  at  two  different  intensities,  3  x  10  and  1  x  10  watts/cm  . 

Both  curves  appear  to  reach  the  same  steady-state  value  of  approximately 
17  2 

3  x  10  electrons/cm  .  The  higher  intensity  beam  causes  the  steady  state 
to  be  reached  at  a  somewhat  earlier  time  as  would  be  expected. 

Figure  20  shows  the  same  type  of  curves  for  laser  intensities  of 
8  8  2 

3  x  10  and  1  x  10  watts/cm  .  The  higher  intensity  again  causes  the  peak 
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Intensity  of  1  x  10  watts/cm  . 


N  (cm' 
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Figure  20.  Electron  Density  vs  Time  for  a  48ji  A1203  Particle 

in  the  Beam  of  a  CO,  Laser  at  Two  Different  Intensities. 


electron  density  to  be  attained  at  an  earlier  time,  however,  the  density  is 
somewhat  lower  than  that  for  the  low  intensity  beam  at  later  times.  This  is 

presumably  due  to  the  faster  diffusion  of  the  electrons  and  the  high  intensity 
beam. 

Figure  21  shows  the  electron  density  for  Al^  with  a  laser  intensity 
of  3  x  10  watts/cm  at  two  diffeient  ambient  pressures.  An  increase  of 

almost  an  order  of  magnitude  in  the  pressure  results  in  an  increased  electron 
density  of  less  than  a  factor  of  2. 

Figure  22  shows  the  electron  density  for  a  48  p  SiO^  particle  at  an 
intensity  of  3  x  10  watts/cm2  for  two  different  values  of  the  absorption  co¬ 
efficient.  Curve  (2)  with  a  higher  absorption  coefficient,  corresponds  to 
figure  17.  There  is  no  significant  difference  in  the  steady  state  electron 
densities  attained  at  the  two  different  absorption  coefficients. 


Electron  Density  vs  Time  for  a  48|x  Particle 

in  the  Field  of  a  CO,  Laser  at  an  Intensity  of 

8  2  c 

3  x  10  watts/cm  at  Different  Ambient  Pressures. 


(1)  Absorption  coefficient  =  3x10 

(2)  Absorption  coefficient  =  1  x  10 


Figure  22.  Electron  Density  vs  Time  for  a  48p  SiOz  Particle  in  the 

Field  of  a  CO?  Laser  at  an  Intensity  of  3  x  10^  watts/cn 
for  Different  Values  of  the  SiC>2  Absorption  Coefficient. 


SECTION  VI 


CONCLUSIONS 

It  was  shown  in  the  theoretical  analysis  of  the  breakdown  of  air  that 
the  decrease  in  power  required  for  breakdown  is  directly  proportional  to  the 
increase  in  the  total  collision  frequency.  As  the  electron  density  is  in¬ 
creased  the  total  collision  frequency  is  increased  through  electron-ion  inter¬ 
action.  After  the  degree  of  ionization  is  increased  beyond  ~  10-^,  the 
electron-ion  interaction  becomes  dominant  and  the  rate  of  energy  absorption 
increases  rapidly,  leading  to  breakdown. 

The  results  of  the  calculations  performed  for  this  report  show  that 
the  atmospheric  impurities  are  a  source  for  a  substantial  density  of  free 
electrons  as  suggested  in  references  1  and  2.  This  is  especially  true  for 

ai2°3. 

The  temperatures  calculated  for  the  electrons  were  relatively  low, 
about  10  “K,  since  the  energy  absorption  was  just  compensated  by  the  loss 
of  electron  energy  to  molecular  vibration.  Thus,  the  electron  densities 
calculated  were  from  the  impurity  particles  and  not  from  the  atmospheric 
^2  or  ®2‘  ^n<*er  these  conditions  one  cannot  draw  any  conclusions  about 
the  relative  effectiveness  of  the  impurity  particles  in  breaking  down  the 
atmosphere. 

The  electron  diffusion  coefficient  used  in  the  calculations  was  too 
large,  being  more  consistent  with  the  value  for  the  free  electron  diffusion 
than  the  ambipolar  diffusion.  This  would  allow  the  electron  density  to  de¬ 
crease  much  too  rapidly  leading  to  the  relatively  low  temperatures  mentioned 
previously.  On  the  otherhahd,  the  use  of  the  Richardson  equation  to  des¬ 
cribe  the  thermionic  emission  for  Al^O^  may  have  given  too  large  an  elec¬ 
tron  flux.  These  two  possible  overestimates  tend  to  compensate  each  other; 
however,  the  large  electron  diffusion  coefficient  may  lead  to  the  more 

serious  discrepancy.  The  electron  density  calculated  for  A1  O  was 

2  3 
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approaching  values  that  could  lead  to  substantial  changes  in  the  total  colli¬ 
sion  frequency,  thus,  this  area  requires  more  study. 

It  was  found  during  the  course  of  this  investigation  that  information 
on  the  material  properties  of  atmospheric  impurities  was  very  difficult  to 
find.  Further  studies  of  air  b»eakdown  should  lead  to  a  great  deal  of  useful 
information  on  experimental  measurements  that  are  needed.  An  example  of 
this  was  indicated  by  the  lack  of  availability  of  information  on  the  type  of 
molecules  formed  by  Al^C^  in  the  vapor  phase. 

The  CELAB  hydrodynamic  code  developed  for  this  contract  should  be 
very  useful  for  future  investigations  of  the  effectiveness  of  other  impurities 
as  a  source  of  free  electrons.  It  was  originally  intended  that  the  present 
investigation  would  include  calculations  of  the  electron  density  resulting 
from  clusters  on  macromolecules  such  as  water  vapor  clusters,  i it  is  ex¬ 
pected  that  the  lower  ionization  potentials  of  such  clusters  might  be  a  major 
source  of  electrons.  Time,  however,  did  not  permit  such  studies  at  the 
present  time.  It  is  suggested  that  this  potential  source  of  electrons  be  in¬ 
vestigated  in  the  future  along  with  other  impurities  existing  in  the  atmos¬ 
phere. 

A  second  approach  to  the  calculation  of  atmospheric  breakdown  was 
started  but  remained  unfinished  as  the  present  investigation  was  terminated. 
In  this  approach  a  Monte  Carlo  calculation  was  used  to  follow  the  time 
dependent  development  of  the  electron  distribution  function  including  the 
effects  of  ionization.  While  the  development  of  the  code  was  not  finished, 
enough  progress  had  been  made  to  justify  its  inclusion  as  an  appendix  in 
this  report  so  that  anyone  wishing  to  follow  this  course  of  investigation 
could  utilize  the  results  obtained  thus  far. 


APPENDIX  I 


DESCRIPTION  OF  CELAB  CODE 

Thi«  appendix  give*  a  description  of  the  difference  method*  used  in 
the  code  CELAB.  CELAB  is  a  one -dimensional  Lagrangian  code  which  allows 
plane,  cylindrical  or  spherical  geometry  through  an  input  parameter 
a=  1»  2»  3.  We  begin  with  a  description  of  the  difference  equation  approxima¬ 
tion  to  the  hydrodynamic  equations.  An  implicit  scheme  is  used  to  solve  the 
energy  equation  including  the  thermal  conductivity  term.  The  method  for 
control  of  the  fundamental  time  step  is  given  and  the  energy  deposition  is 
discussed. 

1.  DIFFERENCING  FOR  CELAB 

The  differential  equations  to  be  solved  are: 

Conservation  of  Momentum: 

au  _  na- 1  dp +q 

»r-  •"  “5m-  <»*« 

Definition  of  Velocity: 

-u 

IT  '  U  (127) 

Conservation  of  Mass: 

V  -  -L  -  i  *R<* 

D  "  a  SM 


(128) 


Conservation  of  Energy: 


Equation  (126)  can  be  differenced  directly  to  give 

J+I  J "2 


(130) 


where 


it”  =  i(t”+1 


(131) 


The  second  equation  gives 

r"+1  =  Rn  +  Un+iitn+i 
j  J  J 


with 


it”+i  =  tn+1  - 1" 


(132) 


(133) 


The  conservation  of  mass  becomes 


.  i  <*rl>  -</> 

j-i  a  M.  i 

J”  2 


or 


D. 


n+i 


J-2 


a  M 


_i  - 


iii 


(Rn+1)a  -  (R?^ )“ 


(134) 


(135) 


For  the  energy  equation,  we  use  an  implicit  scheme  to  solve  for 
the  new  energy  and  pressure,  including  the  thermal  conductivity  terms. 
The  left  side  of  the  energy  equation  is  differenced  directly 
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(136) 


where 


wj-i 


(137) 


To  compute  the  P  ■  dV  work  implicitly,  we  assume  that  the 
pressure  is  a  linear  function  of  the  internal  energy 


P  =  P  +  T  D  E 
H 


(138) 


where  T  is  the  Gruneisen  coefficient  and  Pjj  is  the  Hugoniot  pressure 
which  is  a  function  of  density  only.  For  an  ideal  gas,  P^  =  0  and 
T  =  Y-  1. 

Including  the  artificial  viscosity,  the  P  ■  dV  term  gives 

<pdv'H =  (rDMEH+pH^tpr.it2Qj-l>  -¥*■  <u9> 

where 


<+i  ■  '  v”-i  <140» 

To  avoid  numerical  cancellation  in  computing  this  difference,  we  use 
the  relation 


(141) 


and  the  conservation  of  mass  to  write 


Av"+* 
J"  2 


»♦*.«-* 


2  _ 


<R  2) 

J 


un+" 

-J _ 


(R"+^)a  Un+* 


For  the  laser  deposition  we  write. 


M 


j-i 


At 


n+1 


(142) 


W)"i  -  r 'i 

J  “  2  J~  2 


(143) 


where  F^  i  will  be  computed  using  Bier's  law  and  an  appropriate 
absorption  coefficient. 

For  the  thermal  conductivity  term  of  the  energy  equation,  we 
assume  a  power  law  for  the  temperature  dependence  of  the  conductivity, 

X(T)  =  Xo  TP  (144) 

where  xQ  and  0  depend  on  material  and  may  be  slowly  varying  functions 
of  temperature  and  density.  We  further  assume  that  the  relation  between 
temperature  and  internal  energy  is  linear 


E  -  EQ  =  C  (T  -  TQ) 


(145) 


where  C  is  a  specific  heat  and  E^  is  the  internal  energy  at  a  reference 
temperature  T^.  Let 


so  that 


T  «  £  (E  -  Ej) 


(146) 
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The  conductivity  term 


*c  -  £hRa',&l  (U7) 

can  be  differenced  implicitly, 


The  direct  method  of  determining  the  energy  flow  at  the  zone 
boundary 


is  to  use  the  conservation  of  energy 

*(xiT)g|j+  (149) 

and  the  approximation 


xj±  =  xj±£ 


to  solve  for  the  temperature  at  the  zone  boundary. 


where 


vn+i 

i-i 
J  z 


- - 1 

R°+*  -  RnV 
J  J'1 


(150) 


(151) 


(148) 


78 


This  temperature  can  then  be  used  to  compute  the  energy  flow 
at  the  zone  boundary. 

xaip  .  2XJ-lV|  (Vi;Tj-l) 

l  aRj.  x“7  +  xn*,* 

J  j-i 

Using  equations  (144)  and  (146),  the  implicit  dependence  of  t“+£ 

and  X™  f  on  the  energy  differences  can  be  determined 
J“2 


Tni  -  T?  i  +-57T —  y/ni 

j“  2  j“2  j"2 


(151) 


n+4  n 

xj-i  =  Vi 


1  + 


V*  wn+i 

2c,-4T"-i 


J-2  J“2 


(152) 


Then  to  first  order  in  the  energy  differences,  the  energy  flow  at  the 
zone  boundary  is 


Expanding  to  first  order  in  the  energy  differences  and  collecting 
terms  we  obtain 
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Combining  equations(129),  (136),  (139),  (143)  and  (157)  yield*  an 
expression  for  the  energy  differences  which  can  be  written  in  the  form 


■A  w"  +  2  +  B.  Wn+,2  -  C.  W"+!2  =  D. 
J  J‘3  J  J "2  J  J  +  2  J 


(158) 


where 


A  ,  L  .  Tjrl  TJ-t  .  h-i  YJ-t\  z-  .  -All 

J  '  Y“.j+Yj.|  Tn  4  /  J-‘  cj- 


(1591 


b.  =  it  Jrn?+2  •  av?+i'  +  (-  -1+1 '  TJ~i 

J  J  ■  2  j“2  \  vn  vn 

■l'  VrVi 


t, .  *  . 


_n 

Ti-i 
J  2 


Tn  »  -  Tn 


+  1  + 


Yi  &  "  Y;  i 
J-f  J"! 


l5d\ ,»  1. 

n  I  j-1  M  ,  C.  , 

j-i  /  J  J7  J'5 


(160) 


Ci,Ll^.h^)zn_^_  {m) 


it  5V?+1 


D  =  f"+!2  At"+*  -  (rD?*!1  En  x  +  p"  +  1  +  Pn  J  +  2QnA") 

J  J- 2  \  J  ”2  J-2  H  1  J-2  J-i/  2 

J"  2 

+ (zj  (T^i  ■  Tjn-i»  •  z;.i  (t;- i  •  tu\  ■  (,62> 

J"  2 


When  the  boundary  condition  at  the  origin 


wnt2  =  0 


(163) 
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is  used  in  equation ( 1 58)for  j  =  1,  a  linear  relation  between  W?+2  and 
n  +  ]  .  2 

W,  is  obtained.  This  linear  relation  can  then  be  used  in  equation(  158) 

•  -  1  1  .  1 

for  j  -  2  to  obtain  a  linear  relation  between  Wj  2  and  W  .*■.  This 

a  5/2 

process  can  be  repeated  throughout  the  mesh.  The  general  form  of 
the  linear  relation  can  be  decuded  in  the  following  manner.  Let 


W  I2  =  G  !  Wn7  +  H.  » 
J‘2  J-i  J+2  J-2 


(164) 


Substituting  this  expression  into  equation  (156)and  collecting  terms  gives 

+  C'AJ  Gj-i  Hj-i  ‘  Aj  Hj-i  +  Bj  Hj-i  *  Dj]  =  0 

These  relations  are  the  desired  equations  when  the  equality  is  an  identity, 
i.  e.  ,  when 


G.  1  =  - * - 

Bj  ’  Aj  Gj-| 

D.  +  A.  H.  3 


H  ,  =  - - — iZJL 

Bj  ’  Aj  Gj -I 


(165) 


(166) 


These  recursion  relations  along  with  the  boundary  condition  at  the  origin, 


G  1  =  H_±  =  0  (167) 

2  2 

can  be  used  to  determine  the  Gs  and  Hs  by  an  outward  sweep  through 
the  mesh.  Then  using  the  boundary  condition  at  the  outer  boundary  and 
equation (164),  the  energy  differences  W?+j^  can  be  determined  by  an 
inward  sweep  through  the  mesh.  For  the  outer  boundary  condition,  we  take 
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0 


(168) 


wn+2 
Vi  ' 

which  corresponds  to  no  conductivity  across  the  outer  boundary. 

2.  TIME  STEP  CONTROL 

For  stability  of  tht  hydrodynamic  difference  equations,  we  use 
an  artificial  viscosity  similar  to  that  discussed  in  Richtmeyer  and  Morton, 
and  including  a  linear  term  in  the  velocity.  If  the  material  is  in  the  solid 
state  we  use 


=  -D  [CQ2  1 6U |  +  Cx  C  ]•  6U 
C  is  the  sound  velocity 


and  Cq  and  are  constants  with  typical  values 
CQ  =  2.0  Cx  =  0.  25 


where  6U  =  |^j-6R, 
oK 


(169) 


(170) 


The  artificial  viscosity  in  the  form  of  equation(169)  can  exert  a  tension 
in  a  zone  undergoing  expansion,  6U  >  0.  In  the  gas  state  a  tensile 
artificial  viscosity  is  not  allowed  and  the  modified  form  is  used 


f"D  [Cq  1 6U  |  +Cj  Cb]  •  6U 


6U  <  0 

6U  >  0 


(171) 


83 


Corresponding  to  the  linear  artificial  viscosity,  the  time  step 
for  hydrodynamic  stability  is  determined  by  a  generalized  Courant 


condition 


6R  /  2C0  l6Ul 

4t,  *  -  1+(ci+— g - 

8  \  B 


’  ci* 


2C0  |6U 


(172) 


The  time  step  determined  by  equation  (172)  along  with  the 
artificial  viscosity  assure  the  stability  of  the  hydrodynamic  difference 
equations  and  the  implicit  solution  of  the  energy  equation  guarantees 
the  stability  of  the  thermal  conductivity.  However,  the  accuracy  of  the 
energy  equation  suffers  if  the  time  step  is  so  large  that  the  thermal 
front  is  conducted  across  more  than  one  zone.  The  time  step  for 
accuracy  of  thermal  conduction  is 


A  V. 


(173) 


where  the  velocity  of  the  thermal  front  is  estimated  from  the  similarity 


solution 


_  ST/ Bt 
th  3T/3R 


(174) 


Combining  equations  (173)  &  (174)and  making  use  of  the  thermal 
conductivity  term  of  the  energy  equation  gives 


1  3T 

6t  at 


-L_  =  _J_  J_  fvRa-l  3T| 
At  C6T  3M  [x“  3R J 


(175) 


The  corresponding  difference  equation  is 


3.  INPUT  VARIABLES  AND  FORMATS 

a.  (Problem  description  in  columns  1-60)  Format  (12A5) 

b.  NCYCS,  NMAT,  IALPH,  NEDIT,  LINMX,  INPUT,  NELOP, 
NEPRN,  Format  (8110) 

■ 

* 

NCYCS  =  maximum  number  of  hydro  cycles  to  be  run 
NMAT  =  number  of  material  layers,  up  to  four 
IALPH  =  geometry  factor;  1  for  slab,  2  for  cylindri¬ 
cal  and  3  for  spherical 

NEDIT  =  cycle  interval  at  which  EDITS  is  called 
LINMX  =  maximum  number  of  lines  printed  per  page 
55  is  recommended  for  CDC 
INPUT  =  parameter  controlling  source  of  input  for 
cards.  INPUT  =  0  reads  from  cards, 

INPUT  =  K,  input  from  Kth  record  on 
TAPE8. 

NELOP  =  maximum  number  of  electron-ion  inter¬ 
action  loops  in  ELECT  for  each  hydroloop. 
Usual  value  is  100-200. 


85 


NEPRN 


control  parameter  for  printing  partial  values 
of  electron  loop.  Usual  volume  is  NELOP- 

;.o. 


c.  NT  APE,  NPRIN,  IPRIN,  JINT,  JMIN,  JMAX  Format  (8110) 


d. 


NTAPE 


NPRIN 


IPRIN 


JINT 


JMIN 

JMAX 

TSTC^P,  DTMIN, 


=  Interval  for  taping  time  edits  of  mesh 
NTAPE  =  0  no  time  edits  taped 
NTAPE  =  -1  only  last  time  edit  taped  for 
for  restart 

=  Interval  for  printing  and  plotting  time  edits 
NPRIN  =  0  no  time  edits  printed 
NPRIN  =  -1  only  last  time  edit  printed 
=  control  parameter  for  type  of  time  edit  print 
IPRIN  =  0  print  and  plot 
IPRIN  =  1  print  only 
IPRIN  =  2  plot  only 

=  interval  for  plotting  and  printing  mesh  points 
in  time  edit 

JINT  =  0  active  zones  greater  than  JMIN 
are  plotted 

=  minimum  mesh  point  plotted 
=  maximum  mesh  point  plotted  unles  JINT  =  0. 

TEDIT,  REDIT,  DTEDT,  TMAX,  RTMAX 


Format  (8  El 0.4) 


TST<^P 

DIMIN 

TEDIT 

TEDIT 


=  maximum  problem  time 
=  first  hydro  time  step,  not  used  when  re¬ 
starting  from  tape. 

=  time  of  first  time  edit.  Additional  time 
edits  at  times  given  by 
=  TEDIT  *  (1  +  REDIT)  +  DTEDT 


The  run  is  stopped  if  the  temperature  at  point  RTMAX  exceeds  TMAX. 
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e.  P<J>WW,  WLM,  TEMP,  PRESA  Format  (8  E 10.  4) 

P<^WW  =  Laser  power  in  watts/cm^ 

WLM  =  Laser  wavelength  in  micron 

TEMP  =  Ambient  temperature  in  degrees  Kelvin 

PRESA  =  Ambient  pressure  in  atmospheres 

£.  IZT,  NZ,  THIK,  RAT  Format  (2110,  6E10. 4)(one  for  each 

material  layer) 

IZT  =  control  parameter  for  type  of  zoning 

IZT  =  0,  1  determined  by  AR 
IZT  =  2  determined  by  AV  “•  Ra  ^AR 
NZ  =  number  of  zones  in  layer 

THIK  =  material  thickness  or  radius  in  centimeters 

RAT  =  ratio  for  geometric  series  of  changing  zone 

size. 

g.  MATRL,  MMC^D,  RH<$,  GAMl,  AMWT,  VIC^NE,  SIGIE 
Format  (2A5,  110,  6E10.  4) 

MATRL  =  material  name  in  10  characters 


MMC^D 


RH<^ 

GAMl 

AMWT 

VI<|>NE 

SIGIE 


=  material  model  parameter 

MM(|D  =  -1  Gilmore  tables  for  air 
=  0  ideal  gas 

=  1  Hugonoit  equation  of  state 

=  material  density  at  standard  conditions 

=  ratio  of  specific  heats  of  vapor 

=  atomic  or  molecular  weight 

=  ionization  energy  in  ev 

=  slope  of  electron  impact  ionization  cross 

section  in  units  IT  a^„  ,  / ev 

Bohr 


h.  CHI(j>,  BETA,  ACL,  AEXC^,  VIBC^E,  VIBWE,  VHNUE,  SIGV 
Format  (8E10.  4) 


CHI<|> 

CHI(J) 

\CL, 

AEXO 


VIBOE 

VIBWE 

VANUE 

SIGV 


=  factor  in  computing  temperature  dependent 
thermal  conductivity 
=  CHI(J>(M)  *  T(J)  **  BETA(M) 

X(T)  =  XQ  ^ 

=  linear  absorption  coefficient  at  standard 
density 

=  power  of  density  dependence  in  linear 
absorption  coefficient 

AEXO 

“(p,= “■=(-£) 

=  threshold  for  molecular  vibrational  excita¬ 
tion  by  electron  impact  in  ev 
=  bandwidth  for  molecular  vibrational  excita¬ 
tion  by  electron  impact  in  ev. 

=  molecular  vibrational  quanta  in  ev 
=  cross  section  for  molecular  vibrational  ex¬ 
citation  by  electron  impact  at  center  of  band 


GRUN,  HUGO,  HUGD,  HUGS,  ESUB,  CVS,  TVAP,  WFUNE 
Format  (8E10.  4) 


GRUN 

HUGC 

HUGD 

HUGS 


=  GrUneisen  coefficient  for  solid 
=  C  coefficient  in  Hugonoit  equation  of  state 
=  D  coefficient  in  Hugonoit  equation  of  state 
=  S  coefficient  in  Hugonoit  equation  of  dtate 


2  3  2 

PH(J)  =  C  +  D  +  S  in  dyne/cm 
|i  n 

ESUB  =  sublimation  energy  of  solid  in  ergs/gm 

=  specific  heat  of  solid  in  ergs/gm  degree 


CVS 


TVAP 

WFUNE 


=  vaporization  temperature  in  degrees  Kelvin 
=  work  function  for  thermionic  emission 
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Logical  Flow  Chart  for  CELAB  Program 


0 
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CELAB  (CONT) 


C  SOUND  Subroutine 


ideal 


Is  zone  material  air,  ideal  gas,  or  solid 


Call  AIREST  subroutine 


V  yes 

Are  all  zones  considered  — 


Compute  properties  for  ideal  gas 


Are  all  zones  considered 


Is  zone  expanded? 


Compute  properties  for  expanded  solid 


Compute  properties  for  compressed  solid 


Are  all  zones  considered 


Compute  time  step  for  accuracy 
of  thermal  conduction 


Return 


Competes  positive  and  negative  momentum 


ELECT  Subroutine 


< 


FINED  Subroutine 


HYDRO  Subroutine 


Begin 

__Z3I 

Compute  new  velocities,  coordinates 
densities  and  artificial  viscosity  for 
all  zones . 


Compute  time  step  for 
hydrodynamic  stability 


LASDEP  Sub-outine 


_ L_ 

^^Are  all  zones  considered 


Return 


START  Subroutine 
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CONTAMINANT  EFFECTS  ON  LASER  AIR  8REAK00WN  (CELAB) 


begin  TIME  LOOP  ** 
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APPENDIX  II 

AIR  BREAKDOWN  REDUCTION  BY  TRANSPARENT  PARTICLES 

The  experimentally  observed  fact  that  particles  of  NaCl  cause  a  reduc 
tion  in  the  air  breakdown  threshold  by  an  amount  comparable  to  that  of  much 
more  opaque  materials  is  difficult  to  reconcile  with  a  postulated  mechanism 
that  requires  heating  of  the  particle  as  an  essential  step  since  the  heating 
rate  in  nearly  transparent  NaCl  is  relatively  very  slow.  A  lowering  of  the 
breakdown  threshold  occurs  when  the  population  of  free  electrons  having 
kinetic  energies  above  the  nitrogen  vibrational  excitation  region  is  large 
enough  to  ensure  at  least  one  or  more  of  these  is  in  the  laser  focal  volume  to 
initiate  the  cascade  process.  The  mechanisms  by  which  electrons  can  more 
easily  overcome  the  barrier  to  energy  growth  presented  by  the  nitrogen  vibra 
tional  excitation  region  usually  depend  on  an  air  temperature  rise  induced  by 
the  laser  heated  particle  as  a  result  of  conduction,  shock  compression  or 
radiative  transfer.  Such  a  temperature  rise  increases  the  free  electron  den¬ 
sity  and  with  it  the  probability  an  electron  will  overcome  the  barrier.  In 
addition,  it  makes  it  easier  for  electrons  to  thermalize  over  the  barrier,  and 
shifts  the  balance  of  energy  loss  to  gain  rates  in  the  direction  favoring  easier 
breakdown.  Without  heating,  however,  these  mechanisms  are  inoperative. 

Another  mechanism  which  could  allow  electrons  to  overcome  this 
energy  barrier  without  heating  could  be  the  result  of  the  presence  of  a  rela¬ 
tively  high  density  of  very  low-energy  electrons  in  at  least  a  small  region 
close  to  the  surface  of  the  particle.  These  electrons  might  he  the  result  of 
thermionic  emission  at  the  particle  surface,  perhaps  consisting  of  a  cloud 
extending  outward  no  more  than  the  order  of  a  Debye  length.  Whatever  the 
source  of  these  relatively  numerous  iow-energy  electrons,  their  presence 
could  allow  the  energy  barrier  to  be  overcome  as  a  result  of  electron-elec¬ 
tron  collisions  alone  without  the  necessity  of  air  heating. 


In  order  to  investigate  what  magnitude  of  an  effect  electron-electron 
collisions  might  have  on  easing  the  penetration  of  the  vibrational  excitation 
barrier,  it  was  necessary  to  compute  electron  energy  distribution  curves 
with  and  without  electron-electron  collisions,  and  to  determine  from  these 
whether  or  not  equivalent  densu  as  o>  electrons  with  energies  above  the 
barrier  could  occur  at  lower  laser  powers  when  electron-electron  scattering 
was  taken  into  account.  Although  a  method  for  integrating  the  Boltzman 
equation  directly  taking  into  account  electron-electron  scattering  was  available 
to  us  (Ref.  40),  it  required  considerable  amounts  of  computer  time  per  case 
to  use.  Instead  a  Monte  ^arlo  method  was  chosen  as  a  relatively  direct 
approach  employing  a  minimum  of  simplifying  assumptions. 

It  turned  out  that,  although  this  direct  approach  is  reasonably  effec¬ 
tive  in  treating  the  case  with  no  electron-electron  scattering,  it  is  quite  in¬ 
efficient  for  the  electron-electron  scattering  case  unless  a  number  of  more 
sophisticated  Monte  Carlo  refinements  are  employed.  Resources  available 
to  us  in  this  study  did  not  allow  completion  of  incorporation  of  these  refine¬ 
ments,  and  as  a  consequence  this  covers  only  the  work  which  was  done  along 
with  an  indication  of  the  direction  further  effort  might  take. 

At  the  time  the  Monte  Carlo  approach  was  adapted,  electron  distri¬ 
bution  curves  had  been  computed  by  AFWL  (Ref.  30)  for  the  case  of  no 
electron-electron  scattering  utilizing  a  differencing  technique.  This  technique 
was  subsequently  refined  by  AFWL  to  include  electron-electron  scattering 
effects  (Ref.  16).  It  was  hoped  that  by  making  comparisons  the  importance 
of  some  simplifying  assumptions  made  in  the  AFWL  approach  could  be  eval¬ 
uated  since  they  were  not  made  in  the  Monte  Carlo  approach.  In  particular, 
the  differencing  approach  employed  an  average  power  input  from  the  laser, 
and  did  not  take  into  account  the  fact  that  electrons  oscillate  in  the  laser 

electric  field  only  in  the  field  direction.  The  kinetic  energy  of  a  2  ev  elec- 
9  2 

tron  in  a  10  watt/cm  laser  field  may  actually  range  from  1.  6  to  2.4  ev 
during  a  laser  cycle  when  its  initial  velocity  is  aligned  with  the  field.  Inter¬ 
action  cross  sections  can  vary  significantly  over  this  0.  8  ev  energy  range, 
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and  this  in  turn  could  significantly  influence  the  effectiveness  of  barrier 
penetration,  and  thus  the  shape  of  the  electron  energy  distribution  curves. 

1.  MONTE  CARLO  CALCULATION  OF  ELECTRON  ENERGY 

DISTRIBUTIONS 

In  the  application  of  the  Monte  Carlo  technique,  the  straightforward 
approach  of  following  the  time  history  of  a  single  electron  starting  at  some 
initial  kinetic  energy  was  chosen.  This  electron  was  followed  until  it  was 
lost  by  some  process  such  as  attachment,  or  until  it  caused  an  ionization 
event.  After  an  ionization  event,  one  of  the  resultant  electrons  was  picked 
to  follow  further.  The  total  time  spent  by  the  electron  at  each  energy  (in 
practice  the  time  spent  in  one  of  the  preselected  energy  intervals)  was  record 
ed.  From  this  information,  the  energy  distribution  function  was  calculated. 

The  ergodic  hypothesis  assures  the  validity  of  this  method  of  deter¬ 
mining  energy  distributions  for  steady  state  distributions.  It  states  that  the 
time  average  for  the  properties  of  any  system  of  interest  equal  the  averages 
of  those  same  properties  for  the  members  of  the  ensemble  (Ref.  41).  That 
is,  any  quantity,  A(p,  q)  will  have  a  time  average, 

rT 

A  =  Tin)  •  f  /  A(p(t)’ q(t,)dt  (l77) 

0 

and  an  ensemble  average 

<A>  -J  A  (p,  q)  f(p,  q)  dp  dq  (178) 

where  f(p,  q)  a  the  distribution  function 

p(t)  s  p  (p  ,  q  ,  t),  momentum  coordinates 
o  o 

<l(t)  ■  <1  (PQ,  qQ,  t),  position  coordinates 
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the  ergodic  hypothesis  guarantees  the  ultimate  equality  of  these  two  average 
values,  A  and  <  A  >  at  sufficiently  long  times. 

When  the  distribution  is  changing  in  time,  these  two  average  values 
are  equal  only  if  the  time  variation  is  slow  oomnared  with  the  time  for  the 
distribution  function  to  relax  to  its  q  asi-equilibrium  form.  It  will  be  assumed 
that  this  is  true  for  the  cases  of  interest  here,  although  this  assumption  needs 
further  verification  at  least  in  the  late  stages  of  breakdown,  when  electron- 
electron  collisions,  superelastic  collisions,  etc.  ,  are  becoming  increasingly 
important  at  a  rapid  rate. 

If  the  quantity,  A,  is  taken  to  be  the  electron  kinetic  energy,  E,  and 
if  further  the  integrals  are  replaced  with  summations,  these  results  are 


E  =  “m 
T  -* 


(179) 


and 


<E  >  =  £E.  f  (E.) 

*"*i  i 


(180) 


where  A* ^  is  the  time  spent  by  the  electron  in  the  ith  energy  bin  represented 
by  energy,  E^  The  ergodic  hypothesis  states  that  E  =  <  E  >  for  all  ensembles. 
For  this  to  be  true  for  all  arbitrary  E.  requires  from  equations  (179)  and 
(180) that 


(E.)  =T?  .  (~r~) 


(181) 


Equation  (181)  proves  that  the  results  of  following  one  electron's  history,  as 
described  above,  gives  asymptotically  the  energy  distribution  function  desired. 
Tne  statistical  uncertainties  in  the  distribution  at  any  finite  time  are  equiva¬ 
lent  to  those  inherent  in  the  corresponding  Monte  Carlo  calculation. 
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In  following  the  history  of  a  single  electron,  its  interactions  with 
molecules  were  simulated.  The  reactions  with  air  molecules  considered  were 
excitation  of  vibrational  and  ground  to  excited  electronic  states  as  well  as 
elastic  collisions.  Oxygen  attachment  was  also  included.  In  the  calculations 
done  to  date  only  the  ground  states  were  considered  to  be  populated  signifi¬ 
cantly,  although  no  basic  difficulty  would  be  encountered  in  including  excited 
states.  Emperical  cross-section  data  for  the  various  interaction  channels 
were  obtained  from  the  literature  (Refs.  4,  5,  6,  7,  42  and  43).  The  specific 
interaction  channel  was  chosen  randomly  in  proportion  to  the  corresponding 
cross  section,  and  for  the  channel  chosen  the  electron  kinetic  energy  was 
reduced  appropriately  upon  collision.  A  new  random  direction  for  the  test 
electron  velocity  vector  was  picked  after  the  collision  on  the  assumption  of 
isotropy. 

a.  Interaction  of  Electron  with  Laser  Electric  Field 

The  response  of  the  electron  to  the  oscillating  electric  field  of  the 
laser  beam  was  treated  classically.  By  the  correspondence  principal,  this 
treitment  is  correct  in  the  limit  where  the  electron  kinetic  energy  gained  is 
large  compared  to  a  laser  photon  energy.  The  inverse  bremsstrahlung  absorp¬ 
tion  of  energy  by  the  electron  upon  scattering  from  another  electron  or  the 
scattering  of  an  electron  from  an  ion  charge  center  is  correctly  described 
classically  in  the  many  photon  energy  limit  for  the  same  reason.  It  has  been 
shown  that  when  the  energy  imparted  by  such  a  coulomb  scattering  is  small 
compared  to  the  photon  energy,  a  classical  treatment  of  the  scattering  also 
gives  the  correct  answer,  where  energies  of  only  a  fraction  of  a  single  photon 
energy  are  interpreted  as  probabilities  of  there  being  a  one  photon  energy 
transfer  (Ref.  4).  In  the  intermediate  region  multiple  photon  interactions 
become  important  and  are  hard  to  evaluate.  It  is  assumed  that  the  classical 
treatment  is  adequate  in  this  range  also. 

During  a  single  laser  light  cycle,  the  electrons  oscillate  while  the 
neutral  particles  remain  essentially  stationary.  Ions  can  be  treated  as  station¬ 
ary  because  of  their  relatively  large  mass.  The  temperature  of  interest  are 
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low  enough  so  that  all  particles  except  electrons  can  be  assumed  to  have  neg¬ 
ligible  translational  velocity  at  all  times.  All  electron  energy  distribution 
curves  presented  are  specified  only  after  removal  of  the  kinetic  energy  of 
oscillation  due  to  the  electric  field,  or  effectively  when  the  field  is  turned 
off,  to  avoid  ambiguity.  A  collision  while  the  field  is  on  randomizes  the 
directed  component  of  the  electron  velocity  resulting  from  the  oscillating 
field  and  thus  thermalizes  this  part  of  the  electron  kinetic  energy.  The 
velocity  of  impact  of  an  electron  with  another  particle  varies  during  a  laser 
light  cycle,  and  with  it  the  impact  energy  and  cross  section  for  interaction. 

From  a  classical  point  of  view,  the  electric  vector  of  the  light  in  a 
polarized  laser  beam  can  be  represented  by 

E  =  Eq  cos  (cot)  (182) 

where  t  is  time,  and  the  frequency  is  given  by 

co  =  2irc/\  (183) 

with  c  the  velocity  of  light  and  X  the  light  wavelength.  The  peak  electric  field, 

Eq  ,  in  volts/cm  is  related  to  the  average  power  in  the  laser  beam,  P,  in 
watts /cm^  by 

Eq  =  240«tP  (184) 

A  free  electron  in  this  field  has  an  induced  velocity  component  given  by 


VE  =  V0  8in<wt> 

where 


(185) 


(186) 
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and  e  the  charge  and  the  mass  of  an  electron. 

Combining  equations  (183)  through  (186)  gives  the  numerical  formula 

VQ  =  25.  66  X^p”  (187) 

where  is  in  cm/sec,  \  in  microns,  and  P  in  watts/cm2. 

The  actual  electron  velocity  is  the  vector  sum  of  that  due  to  thermal 
motion  if  the  laser  were  off,  v^  ,  and  that  caused  by  the  field: 

v  =  vQ+vE  (188) 

Neglecting  the  motion  of  all  gas  particles  due  to  either  temperature  or  laser 
field,  the  electron  kinetic  energy  of  impact  with  other  gas  particles  is  given 
simply  by 

W  =  |m0v2  (189) 

where  v  is  obtained  from  equation  (188). 

Using  equations  (187,  188  and  189),  it  can  be  shown  that  the  maximum 
range  of  variation  in  kinetic  energy  of  an  electron  oscillating  in  the  laser 
field  is 

AE  =  0.  6  VPx  W  (190) 

9  2 

where  AE  is  in  electron  volts,  P  is  in  units  of  10  watts/cm  and  W  is  in  ev. 

b.  Interaction  of  the  Electron  with  another  Particle 

In  an  increment  of  time  At,  the  electron  may  interact  >vith  a  gas 
particle  or  another  electron  or  not  at  all.  If  the  time  increment  is  chosen 
small  enough  so  that  during  it  the  laser-induced  velocity  increment  is  negli¬ 
gible,  then  the  electron  kinetic  energy  may  be  taken  as  constant  during  this 
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time  increment.  For  a  given  electron  energy,  the  cross  sections  all-particle 
interaction  channels  are  known.  The  probability  the  electron  will  survive  a 
time,  t,  without  collisions  is  given  by  the  Poisson  probability 

-LvlTj.il. 

p(t)  =  e  11  (191) 

where  the  cross  sections  are  a.  for  a  number/ volume,  n.,  of  interaction  tar- 

1  i 

gets  of  type  i. 

In  the  increment  of  time,  At,  the  probability  of  a  collision  is  from 
equation  (191) 

-A tvEa.n. 

Ap  =  1  -  e  11  .  (192) 

If  the  time  interval  is  chosen  so 

AtvZCTin.<<l  (193) 

then  equation  (192)  can  be  approximated  as 
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is  computed  for  successive  values  of  m  and  for  each  value  q  is  compared 

m 

with  R.  If  for  any  m  it  is  found  that 

q  >R  (197) 

m 

then  the  electron  is  chosen  to  interact  in  channel  m  during  At.  If  after  sum¬ 
ming  over  all  channels  the  inequality  of  equation  (197)  »  not  satisfied,  the 
electron  is  chosen  to  not  interact  in  this  particular  time  interval.  For  this 
method  to  be  valid,  it  is  necessary  that  the  time  step,  At,  meet  condition 
(193). 

The  electron-electron  scattering  cross  section  must  be  included  in 
evaluating  the  condition  (193).  For  this  purpose,  it  will  be  assumed  no  inter¬ 
action  occurs  beyond  the  Debye  shielding  distance  given  by  (Ref.  44  ) 


•h  = 


V 


kT 


=  K 


4ffm  e 


V? 


(198) 


where  n  is  the  volume  density  of  electrons,  and  e  the  electronic  charge.  For 

e  _3 

T  in  deg  kelvin  and  n^  in  cm  ,  K  =  6.  90  gives  h  in  cm.  The  cross  section 
for  electron-electron  interaction  is  thus 


a  =  trh2  =  irK2T/n 
©  © 


(199) 


It  is  noteworthy  that  in  the  summation  of  equations  (193)  and  (194)  the  product 
an.  for  electrons  from  equation  (199)  is 


an  =  ffK  T 
e  e 


(200) 


which  is  independent  of  electron  density.  Since  the  temperature,  T,  will  not 
generally  correspond  to  a  Maxwellian  distribution,  in  its  place  f  times  the 
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average  kinetic  energy  of  the  electron  distribution  was  used. 

If  the  restriction  on-time  interval  imposed  by  equation  (193)  is  much 
more  severe  than  imposed  by  laser  frequency,  then  it  can  be  more  efficient 
computationally  to  relax  thic  condition  by  determining  first  whether  there  is 
any  interaction  at  all  by  using  the  condition  for  there  to  be  one  that  Ap  from 
equation  (194)  exceed  a  random  number,  R.  Only  after  the  existence  of  inter¬ 
action  is  established  is  a  channel  for  it  chosen. 

c.  Electron-Electron  Scattering 

In  the  event  of  an  electron-electron  scattering,  a  target  electron 
must  be  picked  from  the  distribution  for  the  test  electron  to  scatter  from.  The 
second  electron  energy  cannot  be  simply  selected  in  proportion  to  its  frequency 
of  occurrence  in  the  distribution,  because  the  thermal  velocity  is  not  small  ' 
compared  to  that  of  the  test  electron. 

The  probability  per  unit  time  that  the  test  electron  with  velocity,  v^  , 
will  pass  within  a  Debye  length  of  a  target  electron  is  given  by 

2  ir  n 

dt~  ( V  ^2)  =  h /  d9  /  d<P  Bin<pp'(v2>  0,  (p)  \v{  -  v2  |  (201) 

Jo  Jo 

where  n'(v  ,  0,  <p)  is  the  density  of  target  electrons  with  velocity  v  and  speed 
v^  in  a  spherical  coordinate  system  with  polar  angle,  <p  ,  and  azimuth  angle, 

6  •  It  will  be  assumed  that  the  target  electrons  are  isotropically  distributed 
in  velocity  direction,  or 

n1v2,  9,(p)  =  n(v2)/4ff  (202) 

Consequently,  equation  (201)  using  (202)  becomes 
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Vj  +  -  ZvjV^cosp  sin<pd<p  (203) 

o 

which  gives  taking  into  account  the  absolute  value  of  the  velocity  difference 
in  equation  (201) 


n(v2)  ffh 

3 


2 


(204) 


where 


vA=m“(v1,v2) 

vfi  =  min  ( v v2) 

The  quantity  n  (v2)  in  equation  (204)  may  be  determined  from  the  distribution 
function  for  target  electron  energies.  Thus,  for  a  given  test  electron  velocity 
the  modified  distribution  function  implied  by  equation  (204)  may  be  used  to 
correctly  select  a  target  electron  speed,  v_  . 

The  direction  of  this  target  electron  is  assumed  isotropically  distri¬ 
buted.  However,  even  though  it  is  now  determined  that  an  interaction  of  the 
test  electron  will  occur  with  a  target  of  speed,  v  ,  the  probability  the  col- 

M 

lision  will  occur  depends  upon  the  relative  velocity  of  the  two  electrons.  The 
velocity  of  the  target  electron  may  be  randomly  selected  by  first  selecting  a 
random  number  0  <  R  <  1  .  The  cosine  of  the  angle  between  the  two  electron 
velocities  is  then  given  by 


m 

• 

2  2 

L  A(  \ 

3  1  I3  l 

V  TV 

1  2 

V  • /  \l+v2/ 

+  R  |vrv2|  | 

(205) 


The  vector  velocity  of  the  target  electron  is  then 


A 

where  R  is  a  random  unit  vector  and 


y  5 


A 

v1  •  R 


t  1 


(206) 


and  /J 


B 


The  two  electrons  are  now  allowed  to  interact  in  a  center  of  mass 
encounter  and  the  new  test  electron  energy  and  corresponding  velocity  having 
a  random  azimuth  angle  found  and  converted  into  the  laboratory  frame. 

The  velocities  in  the  electron-electron  center  of  mass  system  are 


W1  ~  (V1  '  v2^2  “  w 


w2  =  -w 


(207) 


The  angle  of  scattering  in  the  cm  system  is  found  by  picking  a  random  number 
0  <  R  <  1,  and  finding  the  impact  parameter,  x,  in  terms  of  the  Debye  length, 
h,  from 


and 


2  2 
ffx  =  Rirh 


=  h  yj~R 


(208) 
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The  collision  diameter  is  given  by 


e2/(mQW2)  = 


The  scattering  angle,  6  ,  is  found  from 


6=2  tan 


Equation  (210)  using  equation  (209)  gives 


(209) 


(210) 


cos  0  = 


4h2R  -  b 
4h2R  +  b 


2 

2 


(211) 


in  terms  of  the  random  number,  R.  The  azimuth  angle  is  still  undetermined 

A 

and  can  be  picked  at  random  as  follows.  Let  R  be  a  randomly  picked  unit 

A  -4 

vector.  If  w  is  a  unit  vector  in  the  direction  of  w,  then  the  center  of  mass 

-♦  f 

velocity,  w  ,  after  scattering  is  given  by 


w'  = 


sinOR  + 


Vi  -  (W-  k)1 


(cosG) 


AAA 

sin0(w-  R)w 


V 


1  -  (W  R) 


(212) 


Shifting  back  to  the  laboratory  system,  the  velocity  of  the  test  electron  after 
scattering  is 


<Vv2> 


+  w' 


(2i3: 


From  tuis  velocity  the  electron  kinetic  energy  after  the  electron-electron 
scattering  even  can  be  found  using  equation  (189). 


A  computer  program,  Laser  Ionization  Rates  by  Single  Electron 
Monte  Carlo  (LIRSEM),  was  written  to  carry  out  the  calculations  described 
above.  It  is  described  in  Appendix  III.  It  was  used  to  compute  electron 
energy  distribution  curves  without  electron-electron  scattering.  The  results 
are  shown  for  two  laser  powers  for  which  calculations  were  also  made  by 
AFWL  along  with  the  AFWL  results  for  comparison  in  figures  23  and  24. 
Results  obtained  for  a  number  of  other  laser  powers  are  shown  in  figure  25. 

It  can  be  seen  in  the  comparison  figures  that  the  curves  are  similar 
in  shape  above  a  few  electron  volts  in  energy,  but  differ  considerably  at  low¬ 
er  energies.  There  are  three  plausable  reasons  for  the  low-energy  differ¬ 
ences.  One  is  that  the  energy  bins  used  in  the  Monte  Carlo  (0.  2  ev)  were 
too  coarse  to  provide  adequate  resolution  at  very  low  energies.  Another  is 
that  the  Monte  Carlo  maximum  energy  was  at  20  ev.  Any  electron  exceeding 
this  energy  was  considered  to  have  caused  an  ionization  to  occur.  This  could 
artificially  lower  the  high-energy  tail  of  the  Monte  Carlo  curvet)  and  at  the 
same  time  return  too  many  electrons  to  the  minimum  energy  bin  thus  distort¬ 
ing  the  distribution.  A  third  possible  reason  is  that  the  distributions  really 
are  different  as  a  result  of  the  approximations  made  in  the  AFWL  calcula¬ 
tions  which  are  absent  in  the  Monte  Carlo. 

Monte  Carlo  calculations  were  also  made  including  the  electron-elec¬ 
tron  scattering  feature  of  the  computer  program.  However,  an  error  was 
subsequently  found  in  the  program.  This  resulted  from  assigning  a  random 
direction  to  the  target  electron  after  its  speed  had  been  found  rather  than 
using  equations  (205)  and  (206)  to  determine  its  velocity.  Consequently, 
these  results  are  in  error,  and  are  not  reported  here.  The  program  des¬ 
cribed  in  Appendix  III  has  not  been  corrected,  but  could  be  relatively  easily 
for  this  omission.  However,  another  more  grave  difficulty  was  also  dis¬ 
covered  which  will  take  considerably  more  effort  to  correct.  This  difficulty 
concerns  the  efficiency  of  the  Monte  Carlo  method  when  electron- electron 
scattering  is  included. 
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The  Monte  Carlo  method  employed  in  LIRSEM  is  categorized  as  the 
"crude"  method  in  a  treatise  on  Monte  Carlo  methods  (Ref.  45).  In  this 
method,  the  statistical  uncertainty  assigned  to  the  magnitudes  of  the  times 
spent  in  each  energy  bin  of  the  distribution  are  given  by  the  magnitude  and 
the  time  spent  divided  by  the  sq^'re  t  lot  of  the  number  of  entries  into  the 
bin  when  electron-electron  scattering  is  neglected.  When  electron-electron 
scattering  is  included,  however,  the  fact  that  a  target  electron  was  picked  by 
a  second  random  process  increases  the  uncertainty  by  the  square  root  of  the 
reciprocal  of  the  fraction  of  the  time  spent  in  the  bin.  In  the  first  instance 
acceptable  accuracy  is  obtained  quite  quickly,  but  in  the  second  only  very 
slowly.  Although  time  did  not  permit  incorporating  more  efficient  methods, 
a  number  exist.  Either  a  form  of  importance  sampling  using  a  Maxwellian 
distribution  as  a  reference  function,  or  one  of  the  sampling  or  splitting 
methods  would  be  appropriate  (Ref.  45).  The  rate  of  attaining  a  desired 
accuracy  can  be  increased  by  orders  of  magnitude  when  these  better  methods 
are  employed  which  take  maximum  advantage  of  all  that  is  known  about  the 
answer  ahead  of  time,  and  which  employ  Monte  Carlo  only  were  absolutely 
necessary. 

With  the  corrections,  improvements,  and  refinements  suggested,  it 
is  felt  that  the  Monte  Carlo  approach  to  determining  electron  distribution 
functions  can  be  made  to  be  a  valuable  approach.  As  more  emperical  cross- 
section  data  are  introduced  into  the  analysis  and  as  more  simplifying  assump¬ 
tions  are  discarded,  the  relative  value  of  the  Monte  Carlo  approach  increases 
until  for  sufficiently  complex  phenomenon  it  may  become  the  only  feasible 
approach. 
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APPENDIX  III 


MONTE  CARLO  COMPUTER  PROGRAM 

A  computer  program,  LIRSEM,  was  written  to  carry  out  the  calcula¬ 
tions  described  in  Appendix  II.  A  logical  flow  chart  for  this  program  and  a 
listing  are  included. 

The  program  is  dimensioned  to  allow  for  100  energy  bins  for  the  elec 
trons.  As  run,  the  bins  are  0.  2  ev  wide  ranging  from  0  to  20  ev  energies. 
Cross-sections  are  first  computed  for  each  interaction  channel  limited  to  30 
by  present  dimensioning.  Cross-sections  are  read  into  the  program  as  data 
in  the  form  of  energy  cross  section  pairs  of  points  arbitrarily  spaced  from 
which  mesh  values  are  interpolated.  After  reading  the  cross-section  data, 
an  initial  electron  distribution  is  read  in  the  same  format  after  first  reading 
an  arbitrary  heading  card. 

The  format  for  this  input  data  is  the  following: 

Card  1.  NCHAN,  Number  of  channels,  N,  (£  30)  in  columns  1-5  in 
115  format. 

This  is  followed  by  N  set  of  cards  consisting  of  the  following: 

Card  2.  ICHAN,  Channel  number,  J,  in  1-5  in  115  format 

NJCHAN,  No.  of  pairs  of  data  points  in  columns  6-10, 

115  format 

ITYPE(J),  Type  of  cross  section  according  to  the  code 
(in  columns  11-15  in  115  format) 

1  =  total  cross  section 

2  =  ionization  cross  section 

3  =  vibrational  excitation  cross  section 

4  =  attachment  cross  section 


-3 

PARTCL(J),  Number  density  (cm  )  of  target  particles  in 
this  channel  in  columns  21-30  in  E10.4  format 

ELOSS(J),  Energy  loss  upon  collision  (ev)  when  appropriate 
for  this  type  of  interaction  in  columns  31-40  in  E10.4 
format. 

Card  3  or  more. 

-16  2 

Pairs  of  energy  (ev),  cross  section  (1  x  10  cm  )  data  points 
in  16F5.  2  format. 


After  reading  N  sets  of  cross-section  data,  an  arbitrary  heading  card 
in  80A1  format  is  read  followed  by  one  more  set  of  cards  in  the  Card  2  and  3 
formats.  In  this  instance,  the  first  word  on  the  first  card  serves  only  to 
indicate  the  run  is  to  continue  if  it  is  nonzero.  The  next  cards  in  the  same 
form  as  Card  3  above  represent  energy  (ev),  re'ative  distribution  height 
(arbitrary  units)  data  pairs  to  specify  the  initial  electron  distribution.  The 
normalization  is  discussed  below.  A  run  may  be  terminated  by  use  of  a  blank 
card  for  Card  2. 


Following  this,  control  data  is  read  in  unformatted  input  form 
(Namelist).  The  variables  read  in  are 


WAVE 
POWER  = 
ENGZER  = 
ELECTS  = 

NDT 

NEDIT 
NQNIT 
NSTART  = 


laser  wavelength  (cm) 

.  2 

laser  power  (watts/cm  ) 

initial  energy  of  test  electron  (ev) 

-3 

free  electron  density  (cm  ),  a  constant;  program 
control  constants: 

number  of  time  step  subdivisions  made  in  the 
basic  laser  cycle 

the  number  of  basic  cycles  before  an  edit  is  made 
the  number  of  basic  cycles  before  run  is  terminated 
number  of  basic  cycles  assumed  to  have  transpired 
before  input  electron  distribution  occurs.  (Serves  to 
normalize  this  distribution) 
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NOMOL  =  ignored  if  less  than  or  equal  to  zero.  If  greater 
than  zero  molecular  interactions  bypassed  for 
pure  electron-electron  scattering 
IWRT  =  outputs  computed  cross  sections  for  all  channels 
and  bins  if  greater  than  zero.  If  less  than  one 
this  output  bypassed. 

If  values  for  these  control  constants  are  not  specified  in  the  namelist, 
they  are  assigned  the  following  values: 


NDT  =  0 
NEDIT  =  10,000 
NQUIT  =  100,000 
NSTART  =  100 
NOMOL  =  0 
IWRT  *  0 


Any  number  of  sets  of  initial  electron  distributions  followed  by  name¬ 
list  specifications  may  be  stacked.  A  blank  card  terminates  the  run  normally. 

In  addition  to  the  control  constants,  a  number  of  interval  constants  of 
the  program  are  printed.  Also,  values  of  the  following  variables  are  printed. 


DT 

VELZER 

ANGLE 

TELECT 


basic  time  step  (sec) 

maximum  velocity  induced  by  laser  field  (cm/ sec) 
phase  angle  (radians)  corresponding  to  one  time 
step 

average  electron  temperature  (ev) 


At  each  edit  time,  the  electron  energy  distribution  in  sec/bin  is  given  as  well 
as  the  fraction/bin  and  the  number  of  basic  time  steps/bin  in  tabular  form. 
This  is  followed  by  a  table  of  values  for 


163 


T 


=  time  of  edit  (sec) 

RATION  =  ionization  rate  in  ionizations  per  second  per  cc 

RATERR  =  statistical  uncertainty  .  in  RATION 

NION  =  number  of  ionization  events 

NCYION  =  cycle  number  of  last  ionization 

NCY  =  total  cycles  at  edit 

NEE  =  number  of  electron-electron  collisions 

NEM  =  number  of  electron-molecule  collisions 

NWRITE  =  cycles  to  write 

TEV  =  temperature  of  electrons  (ev) 

EVT  =  kinetic  energy  of  test  electron  (ev) 

VT  =  last  electron  velocity  components  (cm/ sec) 

with  field  on 

VZ  =  last  electron  velocity  components  (cm/ sec) 

if  field  were  off 

NEL  =  number  of  elastic  collision 

NELU  =  number  of  electron-electron  scatterings  for 

which  the  cosine  of  the  scattering  angle  was  less 
than  ELNO(=  .  99999) 

XPBELC  =  log  of  probability  of  an  electron  scattering  in 
one  cycle 

PBELC  =  probability  of  an  electron  scattering  in  one  cycle 
LOGIC  =  value  (T,  F)  of  variable  to  decide  if  electrons 
JE  =  number  of  steps  into  last  laser  cycle 

NATCH  =  number  of  attachments 

ELNO  =  cosine  of  electron-electron  scattering  angle 
above  which  scatterings  were  ignored 
VEZ  =  electron  velocity  (cm/sec)  after  scattering  with 

field  off 


JE  2 


last  channel  chosen 


These  output  values  are  followed  by  a  plot  of  the  electron  energy 
distribution  function  at  the  edit  time.  The  values  of  the  points  plotted  (in 
0.  5  ev  bins)  as  well  as  their  uncertainty  are  listed  next  to  the  plot. 

The  formula  employed  in  the  LIRSEM  program  follow  closely  those 
given  in  Appendix  11.  One  variation  was  a  faster  means  of  advancing  time 
when  there  are  no  electron-electron  collisions  was  devised.  Essentially  it 
first  randomly  selects  a  laser  cycle  in  which  there  is  an  event  using  the 
cross  section  averaged  over  one  cycle  in  this  selection.  This  section  of  the 
program  is  titled  "No-Electron  Channel  Choice.  " 
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pic.il  Flow  Chart  for 


LIRSEM  Computer  Program 
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COMPUTE  BY  ALTERNATE  METHOD 
IF  NO  TARGET  ELECTRONS 
FOR  GREATER  EFFICIENCY 


LIRSEM  Listing 
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